(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.2' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     88285,       1780]
NotebookOptionsPosition[     84061,       1704]
NotebookOutlinePosition[     84456,       1720]
CellTagsIndexPosition[     84413,       1717]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.862841875715396*^9, 3.862841927955557*^9}, {
   3.885758117453233*^9, 3.885758168093186*^9}, {3.8857582077416267`*^9, 
   3.8857582589445057`*^9}, 3.8857584140429907`*^9, 3.885758523608061*^9},
 CellLabel->
  "In[532]:=",ExpressionUUID->"2ba69d38-9fbf-46b8-b9d4-5eaf22c5df35"],

Cell[TextData[{
 "Eqn 6 of Wilson and Head 1994 defines the Gratz number and is (my \
underscores to make the notation work):\nG_z = de^2/Kt= n^2Ed/(Kw_cL)\n\nd is \
average flow thickness\nK is lava thermal diffusivity\nE is volume effusion \
rate from the vent\nw_c is the width of the central channel\n\nI THINK THIS \
IS WRONG BECAUSE IT IS DIFFERENT IN PINKERTON and WILSON, d should be squared \
in first relation:\n\nPinkerton and Wilson 1994\n(again  my underscores)\n\
equation 5:\n\nG_z = ud_e^2/(KL)\n\nu is the mean flow velocity\nK is thermal \
diffusivity of lava, 7x10-7 m^2/s\nL is length of lava flow\n\nequation 6 is \
:\n\nG_z=(nd)^2/Kt=n^2Ed/(Kw_cL)\n\nthey state\nw_c is width of channel\nt is \
time since initiation\nE mean effusion rate\n\nn=d_e/d\n\nwhere\nd_e is \
equivalent diameter of the flow, which is 4 times the cross sectional area \
divided by the wetted perimeter (Knudsen and Katz, 1958)\nand\nd_e<d<2d\nfor \
lava flowing in tubes, d_e = d, otherwise d_e = 2w_c.d/(w_c+d)\n\nequation 7 \
is:\n(,-1-/7)1/3~_.A a 0.o9 (7) where A has the value 0.167 s^ 1/3.\n\n",
 StyleBox["it took me a while to get it, but s^1/3 is seconds to the 1/3 \
power, a unit!",
  FontWeight->"Bold",
  FontSlant->"Italic"],
 "\n",
 StyleBox["E is effusion, g is gravity, alpha is slope angle ",
  FontWeight->"Bold"],
 " \n\[Kappa] is 7e-7, Gratz number is 300, \[Tau] is 2000, nu is viscosity, \
(tau/nu)^2/11 is taken as 0.38 (suggested by wilson and head) \[Rho] is \
density taken as 3200 kg/m^3, 1/1000 factor converts km to m\nNOTE I WORKED \
OUT THE NUMBERS IN THE FIG 8b OF WILSON AND HEAD ARE 10x TOO LARGE! I mention \
this in the resubmit."
}], "Text",
 CellChangeTimes->{{3.8628455820077953`*^9, 3.8628458265759573`*^9}, 
   3.862845915397613*^9, {3.862846018523213*^9, 3.862846104721099*^9}, {
   3.862846138523856*^9, 3.862846339009581*^9}, {3.862846371240184*^9, 
   3.8628464809678307`*^9}, {3.862846542629848*^9, 3.8628466515345488`*^9}, {
   3.8628474062234917`*^9, 3.862847421527569*^9}, {3.862847593224688*^9, 
   3.862847627903718*^9}, {3.8855909656810703`*^9, 3.885590993089445*^9}, 
   3.8855910744556923`*^9, {3.885591127394803*^9, 3.885591332377405*^9}, {
   3.885591720047757*^9, 3.885591752878875*^9}, {3.88572391937379*^9, 
   3.885723930978694*^9}, {3.885763512849751*^9, 
   3.8857635494670277`*^9}},ExpressionUUID->"3f3086f1-6c1a-45b0-9b94-\
de14f6f460d0"],

Cell[BoxData[
 RowBox[{
  RowBox[{"L", "[", 
   RowBox[{"g_", ",", "\[Alpha]_", ",", " ", "E_"}], "]"}], " ", ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "/", "1000"}], ")"}], 
   RowBox[{"(", 
    RowBox[{"1.34", "/", 
     RowBox[{"(", 
      RowBox[{"0.00007", "*", "3"}], ")"}]}], ")"}], "*", "0.38", "*", 
   SuperscriptBox["E", 
    RowBox[{"9", "/", "11"}]], "*", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"2000", "/", 
      RowBox[{"(", 
       RowBox[{"3200", "*", "g"}], ")"}]}], ")"}], 
    RowBox[{"6", "/", "11"}]], " ", "*", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Sin", "[", " ", 
      RowBox[{"\[Alpha]", " ", "Degree"}], "]"}], ")"}], 
    RowBox[{"3", "/", "11"}]]}]}]], "Input",
 CellChangeTimes->{{3.862846936558518*^9, 3.862846941691021*^9}, {
   3.862847038283607*^9, 3.862847043957526*^9}, {3.862847079935561*^9, 
   3.862847080582983*^9}, {3.8628471382885036`*^9, 3.862847168948069*^9}, {
   3.8628472233794117`*^9, 3.8628472953851*^9}, {3.862847724492414*^9, 
   3.8628477494687853`*^9}, {3.885590384133395*^9, 3.885590444424436*^9}, 
   3.8855904883556957`*^9, {3.885590670715776*^9, 3.885590674172257*^9}, 
   3.885590828267941*^9, {3.885590885772645*^9, 3.885590946022293*^9}, {
   3.885591388060526*^9, 3.885591494299769*^9}, {3.885591675859785*^9, 
   3.8855917042367077`*^9}, {3.885591737893363*^9, 3.885591823962961*^9}, {
   3.885591897787376*^9, 3.8855919214894257`*^9}, {3.885591965001871*^9, 
   3.885591966480196*^9}, {3.885641601035142*^9, 3.8856416025309343`*^9}, {
   3.8856417467192993`*^9, 3.885641765005588*^9}, {3.885641919945136*^9, 
   3.885641923262162*^9}, {3.885641959111764*^9, 3.885641984921258*^9}, {
   3.885642682765766*^9, 3.885642688029078*^9}, {3.885643549247025*^9, 
   3.8856435685769997`*^9}, {3.8856441327228327`*^9, 3.885644150178046*^9}, {
   3.885723532428523*^9, 3.885723535286591*^9}, {3.885723839575242*^9, 
   3.885723843502079*^9}, {3.885723897015136*^9, 3.88572391596883*^9}, {
   3.8857525896920958`*^9, 3.885752595138603*^9}, {3.885763507435554*^9, 
   3.8857635106819344`*^9}},
 CellLabel->
  "In[533]:=",ExpressionUUID->"179b6fa4-09df-4050-896c-1e873474453a"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.885591971174508*^9, 3.885592004730097*^9}, {
   3.885641880763339*^9, 3.885641887889543*^9}, 3.885642009572893*^9},
 CellLabel->
  "In[534]:=",ExpressionUUID->"c4ddf92d-7740-4248-9d94-e0085be039d5"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"L", "[", 
    RowBox[{"3.72", ",", "0.1", ",", "\[Theta]"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "1", ",", "10000"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellLabel->
  "In[535]:=",ExpressionUUID->"939136d3-589b-4ee0-8a45-d858bfc87826"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVlGs8lIkChyfLD9mSTjvk0hDlth05ytaJ5q9xO0cpU+M2G40pmTFmeHMi
kdhyxkbFQbUhurjEIiZDimZeueQuW0cuS+XS7OlGMYi28+H5PZ+fL49JkIh5
RIVCoQR95f/+Nc157jfKezo1RLzTVN5IX3It2xbjvApP1Iwnf/cboV/kaU3X
iw0hNvq7rfXYa7qOq+RpQrAZ+swHS8z9PtD9n72wuym2wkv15ZXtiZ/oTInY
ncvajEvhNgUWp+foM59pUZHBdhjm1wWrsRfp9POefVni7aga6NWJP6SCV0/H
VIquOKC1/F7hrKUqxLQ429rbdGQbxcZkKdTA+y+rOWx8FyjjjElDliYqeWfM
Lmr8A21HJzpqRrWhurK7Ti/tn3A/UZG36/4qeFca7M9fuxsx/5mfEZ/VwcJ8
ZUKl1V54X2hqdVz9F6xvChfxC5jwnW7rdtKlItTnWntKiy+aXxx+GzFigAxL
F72eFD8kHe8SpWUYon5Bwf3Oyx9Gr9u+qXYzwqo8+4Wcfja+SbJhbS5ZB8lk
h2W5IgARmj7KoSpjfHE4RYYkBaH0tpUHI8AU7p2lxueKQhCmdMwhTljgyWoT
/tvVPOyN3LOuussCAT6ZVV5xPJza9tJausESkSNxrnpMPrrOxi7NdFsib8pT
cGshFMOWrdnnzayhpH6QNniK8HqmbF5dtgl95b+Yx8wTcOvgVcfzbXFWOVxf
lncCp4/OWGc52+PFmF5mvdpPULXdwtNu3IlfrkvdlAbJmOK033XyZKCWKtji
XH0B0AhPjwpww6vf5IGj79KxX3l1fErNA7RYf3mxaxYMtK1ls4OeGG8TcAdS
r+Bzu76KobMXNpTe2V+QnY2mP+T6/u+Z0A1Zxol5cw3RA8GCUOkBOD2XDU8n
5iMyws9H7OsNRxa+D1S/gZhnlpNKhQ9W73TW6Cq9CVpNKmPNFT+sF42e07Iv
AMXUbL7LjI1vV6rm7h4pxA4/8fK/tf+I+8sO75UnFkOqI5m/ERyAIseDRaqa
JYh03xdlMhoI4cTABY+KUiRtHD94s+4Qwvun1wXuKINAoynQNYoDj+1pDlET
5fDlOHGuIwizM3c/jz6qwI3xkeYoVS6CP70re5t8Bx9twvZcruVCfltQFBdQ
iQ4aLZMdexgOqtWfpsyq8Oia/b/km49ATyM3PXGiCiskvWzWyBGMyluOZ9yT
oLzb/LEwJxgtJg/LUuLv4nzJ/m0lu48i+4f5TeqsahhTqN57NUPwZLb389WN
UhgeY3czZCGw+DGapquQgnqq43pzKA/8qEPJa+pr4Ja/LzfDlI/H9t+5KX+q
xZjw4fKNnXyk5oUcmHG5h/ZK+VUiMhSOsqWi43p16O4SiBTmAmglcbyGBuug
cnrufHSnAE3L1J2a7txHRH7oIWpiGK6ZnDiTGf0As8NmmR8thBj4I1k3fkc9
ahuHDZvbhFANG/qgWNkALit+NuWkCB7pizTHvgYM3rZU8LXDQf+h4rEG6yGS
l/O4sUfDcW5+8yZK/UMMuITapD4Ix7j25VoFTYbWoZ974lZEQHo6j0f9WQbd
iYWFDm4ELjEuMLa+kcHOirtrhSQCxew9LhxPOXK0WPF8FQK2XaX8TqkcvpHm
Kfo+BFoypmr/qk+C86j/195bBHZrvV8aMyBxa4P0jG0hge6EN4wcIxJz9s9Z
F4sI9AsnO7RMSEw+nx7fU0JA4T48OmlOotw0qbGlgsCKxVbNG1tJOG8pC3tQ
R4B5JN9P14tEwxen2Os9BJ4O5OZ2MklEXabkUJ4Q8Gdmvzx7gETYWr2agD4C
XPqlsI8+JGzsDgzpPyMQuTY1oSeABNtP51X6IIGsjujicwISuiWmLxMmCOg7
H3+7S0hiX4W0+PdJArn3jtnNi0jMfFvGd1QQKCwUPgg5RiLIrHFA+T8CNQmH
e1xiSJh8XxAjmiLgoORQF0+S0PQ8uKZjmoBMGMiWxJEwiDxZZPWJwGO2/9j6
hK+9xu8axmYJePb6WD1PJHGqx4XBmCPQ684SpZ0hMdavJsv7+g3vBqbEPYmE
rNNm+9ICgYGt++a+/JtEjnNjCXuRwJ/9v+/5
       "]]},
     Annotation[#, "Charting`Private`Tag$50703#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.885642023012013*^9, 3.885643498452956*^9, 3.885644159171473*^9, 
   3.885678891558364*^9, 3.88567908745133*^9, 3.885723549118505*^9, 
   3.885723640950344*^9, 3.885723730588478*^9, 3.885723936121622*^9, 
   3.885751510036446*^9, 3.8857526104307127`*^9, 3.885763563197587*^9, {
   3.88576676848603*^9, 3.885766791258297*^9}, 3.8857670585991573`*^9, 
   3.8862454684143047`*^9, 3.886245585912067*^9, 3.886245656724668*^9, 
   3.886245724728736*^9, 3.8862458099226027`*^9, 3.886266173229937*^9, 
   3.8862662057607527`*^9, 3.886266831422563*^9, 3.886365643028349*^9, 
   3.886366037971958*^9, 3.8863660734047823`*^9, 3.886366129747382*^9, 
   3.886366228882841*^9, {3.886366264610959*^9, 3.886366347592742*^9}, 
   3.886366593155428*^9, 3.8863685845748796`*^9, 3.8863686925986967`*^9, {
   3.886368816717477*^9, 3.886368830873003*^9}, 3.88636894644285*^9, 
   3.886368980650627*^9, {3.8863696615449047`*^9, 3.8863696861738358`*^9}, 
   3.886372266100079*^9},
 CellLabel->
  "Out[535]=",ExpressionUUID->"55a4c815-d58f-4dda-9775-7cb5ab561b3e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"L", "[", 
  RowBox[{"3.72", ",", "0.1", ",", "1000"}], "]"}]], "Input",
 CellChangeTimes->{{3.885641386939507*^9, 3.885641399756242*^9}},
 CellLabel->
  "In[536]:=",ExpressionUUID->"df530b00-027d-4520-bd6d-95f51302ed02"],

Cell[BoxData["46.18020101314635`"], "Output",
 CellChangeTimes->{
  3.88559200551821*^9, 3.8856414119427013`*^9, 3.8856417799950943`*^9, 
   3.885641895741084*^9, 3.8856419702851686`*^9, 3.8856420230218973`*^9, 
   3.885643498478518*^9, 3.885644159188431*^9, 3.885678891565886*^9, 
   3.885679087459178*^9, 3.885723549128031*^9, 3.88572364096023*^9, 
   3.885723730598481*^9, 3.885723936148321*^9, 3.885751510047089*^9, 
   3.885752610442904*^9, 3.8857635632222567`*^9, {3.885766768495853*^9, 
   3.885766791268539*^9}, 3.885767058607689*^9, 3.886245468429535*^9, 
   3.886245585921269*^9, 3.886245656733616*^9, 3.8862457247380238`*^9, 
   3.886245809932638*^9, 3.8862661732705927`*^9, 3.886266205774568*^9, 
   3.886266831434689*^9, 3.8863656430416937`*^9, 3.88636603798342*^9, 
   3.8863660734141207`*^9, 3.886366129760665*^9, 3.886366228893161*^9, {
   3.886366264619981*^9, 3.8863663476020107`*^9}, 3.886366593165181*^9, 
   3.8863685845849733`*^9, 3.886368692607576*^9, {3.8863688167260647`*^9, 
   3.88636883088043*^9}, 3.886368946450965*^9, 3.886368980659492*^9, {
   3.886369661555126*^9, 3.8863696861843863`*^9}, 3.8863722661097403`*^9},
 CellLabel->
  "Out[536]=",ExpressionUUID->"9c087279-646f-4c5f-b2f6-56516af28d56"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"L", "[", 
      RowBox[{"9.81", ",", "0.5", ",", "\[Theta]"}], "]"}], ",", 
     RowBox[{"L", "[", 
      RowBox[{"3.72", ",", "0.5", ",", "\[Theta]"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "1", ",", "10000"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.885677410699976*^9, {3.8856774750991898`*^9, 3.885677515957527*^9}, {
   3.885678847267879*^9, 3.885678850796009*^9}, {3.885723596420101*^9, 
   3.885723603523842*^9}},
 NumberMarks->False,
 CellLabel->
  "In[537]:=",ExpressionUUID->"7bc97526-b8d6-424c-b9e7-d44e6eff27c2"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVkXk4lAkAh9XahnVbIbVZEk2ltCieML/WEUtkHOMIOTPjmPmqVdmdFSnj
6BhRqsnRZUwi56Qpx3ztpHbZVMpRNlmUXWyECWnbP97n/fd9ntcogk2PXqyg
oBDxmf99g+/0oVPhX9rbicD2LrmUtuBSYZPspAl3IedwwMBL2immymQjbwUk
z3UuikXDNC2X2mepMSaoGB0m5XXjtKDnry2v8NZi5FyOxphwkkav5blG+llg
5rS1Q4znDG163vDA/hhLDL6foTKeztJoJzyfnuHZ4mFWEV/0mwL+eja4WHjO
Dhp1jcxfExaDZ8jd1CCiYUdR04EHFEUwu/zuJwx9j51XX91TN6OgmpluckrJ
DUf4lXZRq9SgqP5Ios//ATe0amw2nlWHf/Vyn5JlHlgVsEsx5wtNzM1Wp1av
9cK7+xo13yRqwVjGYbOu0VHeaSMTWOggjlH0e05rAITjWziM5GXIozrrd+QE
4ojeIvn0SgM0zo1ELvUOQsd6g7VVUgNoFm+eu9gdDK8s89N5lBWofdNGrRwJ
haTzQvXzjJX4ZPcLGXssAh2dFmkFNUZwbS//NlsYi4JbXTz3UFM80TZijWkz
cf30Qr3xDVOEMvJrvLlMDIp41Ltzptj/iuuiT2dhSd0UFwVmKJ7wjL86F4f7
jxxNi9rXQK77TtzkyYYLUjJzN63D08rzZsmzeyFxSa0p79yIo/K+xoriQ9h9
2OLEV0FWeD2on9/45RG8VHdp8Vpph/OXxNvlyzOxUcDNmLmyDQ268VZO9SdB
C+QxRMHO+KtTGtY/ngsf+mD91WpXGP4cJC1zOQN6gdWer609MPRbfGTv8XOw
nL1ZfvikF1aXV/lcEwiQ3NaWa+LlDb3YReHJo0UYph7VVxunY1tPS99kWgmW
fXfU1zneF/Z+WB9GuYwXzvP7dF74QdvBSemP8iv4xfTn/vXRDBiz+7NVNl9D
oaByoKc9AKrqioUer0rxPk5u0v75w51FUV7StDJcWC1yG2sNhtA+RKiofB1K
xPTf3cYhSBzuPel+sxxWfT3jamdCwemeXBm2tQK1lzqMKKNhcLfl2x0YrkQU
654su2I3Zqbr5vt/vQlxKbEiIzwcMVPjFWOZVWDt5WZpUSIgFcULuaHVyG5+
y9BuiICdYv3UhEkNVpWaGkyGREJfqTA3bbgG5+e4DqyFSPRLW5PybtfC/lXJ
lP31KLQaNVfkpNThyQJTO8Y9GoIts+YUv3rsXuzPb5+JxpOZx/MXTMWQLB0u
HSuKwZpdBw31RsTwbqikyrbuAevA7kydxlsY68sykw7swcPNS7fLjzSAeqlS
71pKLI4Xx/pOO9/GQTWbLZf1mLBvWRAm6UsgKNZNKZMwoXIs3PvlCwkeyBZ4
r31YkC2ibJNV3cH75oZMFTkLRUaH0vMP3oWbztV7x0/FoffvTL2UrY3YWXWL
CFkdD8WEl+9G1JuwTlho9oUsHu65Hw3tnzbB3nbJnENwAmhbbj5U8muGuY+i
zHY+AdmzFuYKjc0w0HDfm30uEUMaBQ0jhi2IHhDfEK1jQ3y4mKmb1QJa0D8b
OA/YOOt40tF6tAWHGMr5Ghs4KAve4RzuKUU6RTO94kcONv1RzmoXS3HqhHRK
h+SgNW+iYYMBiT+FNjNLVAh4qPy7MLicxJq4XRe3qxJ4lDrqePEbEuaqTg48
NQLdiW/aVIxITBiHc5Q1CYy49vW/MSPxLfV1vspSAmofHyhftibxmPXTkLoh
AXp0SaCeNwmnuoESLUsCz3oLC9vpJAqFeSXeVgSC6IKBo74kTn+1XMC3JhBJ
O5vwnkGC6ixK1bIhsH/Z8dSOUBKbJf3qWvYEzrQdLMuOJ1Gabv6dxnYCBk5J
Y98nkrC1SPrT05VA4e19lrNsEs1vBRkn3AiUlibejd1HYkh0rFXNg8Ct1KgO
52QSdiE9qqreBOzk4boffyIxyDHmu9MJtCSGBddySTh2eWhm+xB4GBw0aJxK
4kvyh0/K/gQ8HzPW9qR97vc3ZLsxCDx29WPz00n4hnZ18wII+DfRa12PkVCI
SKK1BhLotd754VMGiYi3H4oowQT+A2rZ4EY=
       "]]},
     Annotation[#, "Charting`Private`Tag$50751#1"]& ], 
    TagBox[
     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVk3k81AkDh7VpqUFyzqxtZ6c3hdLhSjnmaxvH26E1GWqUc8m4RgyD5rch
wqLDql2FlIh4J2nQKMfMb5PaEKXdUt5UyrIhJuTc9o/n8/zz/Psw/PnswC+U
lJT8P/Ov/3ea9albaZRJUdENc6A9ZM47ia0TWJq41pfXdezAX8xTPMp4Y/rX
KPX3f03VVTBXOEmeJAWtxuyKgn7hwDST+8cr88vpJhCaLT1eIVYCW5LuEsDZ
BBONgX55y2JMzNKFgiBzSJdXRtJ++xLME66Pz6ZvRRV7ouXvRDW8edL/RVme
LeZ5BSUip+VIpxObpVeZUGgE0YwXNMH7k3M3/O13aFjP/t23RBvVvJTVp1T/
C7k4S+uSlApljYe3qKd3oNaee2flRxo8qg32XqTtQl11Zu9uUwPMTFcnVZvs
wUrHhaGC3JVY1RLJDyllwzCYntzGYCDU88KDrNZ9CBcxBxILDZFr7EjtzNqP
/eMb8jebrEHjzGCArhsX4TTnWnrNGmgWWc0UPPWCH4US2d66FpKBNuNrg95Y
pnlOjzZgjAXbH8ng4/4I7T4gS1u2AS7tld9mlgWD0paQt2TGDI+0GCHDWjzc
c7UZG7Izh7fnmRtuBA9vQhjywkRzCF4STlR2CKrevq4JX2KBojHXsJKZUHR6
GZ67r2aJKb0PdU2ufBQMfZqz1N6Cx9fOrU2YjkLJpFOx92IbpE71NoqL4lG1
b2FewXLAq37qmcYlx+Ac2u1Efe6Mc5fqnKcMMhCf7LJR8mE3pHphFqzak2hv
8NE7RWfjTbfcp28kB1buwgfqH91BF3Hl5U5nkVjD2M0b8cTb38MCerLzkNet
tu/RXS4MK6/vLc3PxzKmTby670HoBy/yS3h/AbIInZQPmr5weCbrHU++CIOD
xGSYux/sOFjvo1IMsmc9Ky3LH1r2LNWOyssQmF2K0XsWgFX8vkyKVSnatu3Z
OacWCDUN5cJdL6+gcJtxLocThNuLftgjTy4Hj7hqXfLzIZTZHSxTXloBf51i
yXppMCLe9ZzcWVUJ043Nu4wUPEQ+Hf/Gx0YMZrjs/Eb9UOzcetpW+O4a/qgY
/T5pTxgmJ2pm++5UQVs9oh9HwxH0cUQ8nHEdEeL8sU2/RkB+NayM8K7GiDGF
69zCh61y7cex1TeQlcStp/AjQVUtzEl+dwPaaVVMy8lI9MlbY3PrJejKfqH+
POMwWhnN4qyjNfgzve2YfEUU8rdMm6pwalFbPeihcTkKjya7Zs+vqUNv53RH
nVk0jA7E0fUH63A7XWtSURuNEKFvhk7jTex9NNV+1FGA+1a6zlPHpDDVPZx8
954A2UXB7hOO9RhjSyqadsTATjZfFku9hfPvczjD3TGgHPdze/H8Fm5SDD1o
nrFoWaTi0HL9NtrmC2IDX8XiAiM+5UxcA4itgTntAUL0DGXoH7VpxEJjjEjn
lRDK4S8+DGo0IfDnu6mBIXHYmTNHt3vchB0T4xYqf8WBuaXqviqnGQdLUiOL
Dscjc3qTqVJjM+503zzyzWQ83i7/VTpIlyF4Tb3119EJqEss4un9JEOc/b3s
2akE/LL95HbL9zJMrlJ1FMQcQbnXbkc/VzlSpSpspZEj2NxRGdJeJ0fVf9bZ
ekWJ0Jo7Jt3wFYlvdexZgr9F2EUZne83IFEz/ZqfOyzCw6T32wtWkrh14kqO
ZFSEpxEDbRQGiUUpmQ0KhQiDLr19A2tJaM+8LBPMiaA+d29psSUJI2bukygN
AuzAi/v13Ug4lRfnBm8i8KSnsLCdTaJwdUdJmhkBLjv/dao7CVcDK3GpBYEA
5i/hCk8SQUb1l99YExDQspM6vUlQhBIjXwcCZ9viyjPDSMTN22dz2AS+YsUO
fxdBgr7ulJbAnUBhfbT5NJ+EZcy6kzkeBK5ciWgIjiahpsyK6OASuJn0Q6dj
wuc+Wu83lwACtlN+enNHSFiMdy4+FEh8/szHS0KQGHrYuy31EIH7Xtz+VUkk
FA6jJ2ShBFy7PE2eJZNA8UTF/8MJdLlw+KdTSGjNuzXP8Ql4NLElLsdJTFkv
fmAQRaDH8vtPC2kkPLZROrYKCPwDHy3m3A==
       "]]},
     Annotation[#, "Charting`Private`Tag$50751#2"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.8856415795761547`*^9, 3.885641780015531*^9, 3.8856418957582006`*^9, 
   3.885641970303458*^9, 3.885642023039548*^9, 3.885643498488956*^9, 
   3.885644159208355*^9, {3.885677505724729*^9, 3.8856775177215223`*^9}, 
   3.885678891597692*^9, 3.885679087488453*^9, 3.8857235491604424`*^9, 
   3.885723640998955*^9, 3.88572373063844*^9, 3.885723936191918*^9, 
   3.885751510082061*^9, 3.885752610509244*^9, 3.8857635632549458`*^9, {
   3.885766768533608*^9, 3.885766791308795*^9}, 3.8857670586494207`*^9, 
   3.886245468465057*^9, 3.886245585956362*^9, 3.8862456567672663`*^9, 
   3.886245724768363*^9, 3.8862458099695663`*^9, 3.886266173325808*^9, 
   3.886266205814196*^9, 3.886266831473998*^9, 3.886365643077194*^9, 
   3.886366038018132*^9, 3.8863660734472237`*^9, 3.886366129813059*^9, 
   3.886366228929805*^9, {3.886366264651676*^9, 3.886366347640233*^9}, 
   3.8863665932038527`*^9, 3.886368584621955*^9, 3.886368692641343*^9, {
   3.886368816764702*^9, 3.886368830911436*^9}, 3.8863689464847527`*^9, 
   3.886368980703849*^9, {3.886369661599715*^9, 3.886369686229147*^9}, 
   3.8863722661561117`*^9},
 CellLabel->
  "Out[537]=",ExpressionUUID->"a66dcbea-304f-4b4e-98b4-ebe2997a1f88"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"L", "[", 
      RowBox[{"9.81", ",", "2", ",", "\[Theta]"}], "]"}], ",", 
     RowBox[{"L", "[", 
      RowBox[{"3.72", ",", "2", ",", "\[Theta]"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "1", ",", "10000"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.885678867056327*^9, 3.885678871291081*^9}, {
  3.885723610627496*^9, 3.885723619941792*^9}},
 CellLabel->
  "In[538]:=",ExpressionUUID->"28d0548e-1f8d-405a-9184-5bb1b9d5240d"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVyXk81AkDx3F6ppeR3bTK0czKkk08KTZHMprvxjhWD5lIEXKucYxrMD8t
ratGpO3AVnK169hlXJMJzYjfFl2DcjylDWUku7XLbK6k3f3j83r/8dEPimGH
rlJSUgr6p3+tO+uwOKj0J/NW3NNEWsM95oqjcFeKwzoEnZTpWPTKmd9x1BRS
wacwu7RKT501w/zEUTSUHmaIB9NfFDVvnmf6DD/f+YPABEcKi8I2NC4z2SKB
c7CXGUT1dfcblZUxt6yXzAvbifWHt/QMn/8PmPluA4UCG3QearioyKBiYki+
qvoiA/WWiy0G8jUQ6KWat/7ERL9dx7VM24/B+b9Xd/TkXoT3nypdnliHJk6W
4XdUF9hSXx80sNcEZW1fu87Zr+ASk1jJZWvhYBP9QPnGfZBef5sh99fGu6Wm
9CYTdxjqjqpSYzfC4HZsTEQlG3fWHgjRSdZFpHfp/byeQ9hhYfTcr8gAF4xZ
Ov15h6Gc1kdPU98M6bvpYE0PHwRLXyzPndiMdWVW76489oX7S8EhVZ4hRFMP
jOun/VE2ympjOG3BB0YaGX4iCCnHWDu4g8ZwltV+llsdDqHHYHgq2wyPNPQj
3mhwwA2u2bQ6zwz+3gXNHqkcnB/2yzx6ywy8sVRHHXYEQgPethVYm6Ns1i3q
x3eR2PrKbF6d9gUWtGbEHW4x6J1WHRkd3omB+ktGKUvx+MotTTa42xrZC8+k
wjICFY6/ZOY32+G5XKdAujoTbbusrpRmOeBShdhpgZ6D2WsT0g03XdCqFWXh
0HIG+d3OQ4L/umNisCtg/I9zmLBkonI1G3rf+HTVOBYiI45UZxR4YvJeVPDI
6Yto85NHlKt44/PaxgOVxcWIuzrKGNh7GNrhyoEpr0thtkw0ycp98eWTzmeK
jHKsOTlkYyzzg50XtgWoXIX1p8u8mckAaOxxoPbW/gA5RRS9XyMQBjHjuWpW
lej0H8y3pQXho7WUkn1jVfBUEO5FFsG4oRzi3pVRg9RK121sxxBU2/lVU1R/
ht2PNI5kbyi4L0fOuDbUQvKKcuayRxhiHys2BdgK8erOI6cx16/hanOWkfyy
HhVTzmq8Q+GYn7u2PH6rAQU2tUvVfhyEvf1D+CanES5PzffAPQJdP0VVp/o3
gXJkJE/lQCQYlJa3s4bNKFZaNK+yj4IOteRcxstmuFp0s4ZMozHe1ZN0oU0E
66OLEnUTLnr0bwrzjl/D1/TOY/fWx6DYeslUxasFGnEKhkglFo/mHy5f3iLG
1pQ5YVBGLLYe4etpT4vxyZbwbNpfsYhIPpqzQXodAx0JZUJuHO5aaTotZLZi
cIdGTsxYHE6XhXvOsdrwnsK+7O0eD7vOleoknXYMa/orU7vjoXYi0OPXp+2Q
SdTv060ScFtZ5cvbjTdgEeopkNcloFSfyCrgS3DfNal9uy4PI7/laB+3laKv
XrL7Vg4PlOhfZ6bXdiCLM9CxssKD67n3enYDHTjVPW48EZkIpnXDXarXTST/
LPzFdzwRuUtmpkrSm4gVK04beSZhUv371mm9TmTSPRq7pUkQf1vG0TrViYTs
EKUp82QU2Z+xt3zdidGi76kzxcmo8f0fK9CtC62PEnWTV/Nh3lsbIRN3YVPC
ybCoJD56Lsy2bqeR4GVfF0eO8rFP7c8VOZ1EUnphLXOcj7701/ZXdEk8aeq7
uP4FH4+5Uw/U9EnkKo74tU/yMe38bHzKiET7+d8Ea97w8fH7O6pXLUm0Ugrj
q5b5YIeWH9b2+OdrSwpGNhIYGikpkbFJjNU1FAnpBHzYxS+yPUmIaLRz6boE
gplF0X95kzCt8Ig10ifA23g6vd+fBDOL8zB+K4HCB/ya3CgS9YrfTanWBGgO
SW/2cknQc4rFT3YRKGlL2LkUQ6LilMymbjeBqiquJDyBhObsDSP2HgLX00P6
WSkkjps8rrnMIsBYCNR6f4zEwGLjfLQTgU5ugK8olURa3Po9cCFw19dHbpBO
QiUxQTSxj4DbQ2+TJxkktuUnTLS4EXjo7BVzNovE5xEfPsrZT+BgB1vkfIKE
da/2dl82gRHL/YsfTpL4NlPiZOpJ4G+ibuQp
       "]]},
     Annotation[#, "Charting`Private`Tag$50800#1"]& ], 
    TagBox[
     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVxX880wkDB3Cih9KpPGJFuVGPm7tSKelp2sevrIispLgIha0xIcMYw1c0
PcVZUfKzi0pE8rOm7SscceWcTtKlSEdKU25Inuf54/16U/14rOMLVFRU/P7n
/9/KsJ/+XWWC4biYoZYlfMr4uqvcKsZ+GcQJgu8uiRSM82ytSWmqIWydLidN
uMwxlu+q7hUFrMW4vpr0dJEaPJ++sriaaoYvQpujRYc0wapOZfq7b0TCx0NF
V54vwdQXI35EgAWsFl7pbcldBsZ/XHoupG6H4ZPVU1mPdDHUO7ygNIeOn3nT
nYe79ZBqFLep/gYDQt5X61VdFLD/cG8NfmMLU7KDPXzTEFXs5LXnNXdDPEy3
alpsDHXtx42UjD24fHLnjU8TxjhYZbC/cKUzuFXZqfIeE8zOVImqzFyxNPiq
jUyyDsYtoTzONRac5Zt7S1VoOOGR/yi97RB6rd4t1z1rjiyaA+VJ+mHsNhFU
Rf2wEdLZUf8Vbp5I4D/26svYiGUFlrNX+rxwy+LUazWvTah+20mrGPVGDy3/
UcTIZszThWRQih+0kgrfdYxtBbOr7FtxaRAS23WXON7Zgd90qJz3Omxo0ptd
/9Klw9tDcsctjo3tddIjiKQj4mXcLgqLgx2NCRIdK2sUKFy4P8+ewAOf+nfc
up1Q6n2sbXLhYcD44PeiXqCn4pJpzEwY+A7aarQaOxDKF9LygmhovDnWMGWx
G6+GKRLpwiTcA2uxPtsNl4pqHZUGaUjXNkd5nzvq9bhb7GvOQas5SDQf5omh
3+U+gx8y8fm5B7tN4g2jWE/59V0XcHtLDM+c74s3HVz//rM5+IcsXXH3jD/W
lVXuv5abi8BLF7UDeMehH6TqGzOeD9sUYVpYWCBsnsleTCYWYt80bTxtDxvW
7vjBR6MYe99f61uocwI6O+01fy27ipjQ0G+7lVwY8wbFWpbXEE8pGampC8ES
bfU855cliI9zFE06heKe6jFXeeJ1VCibLbdvOYlS6yOl6otu4vSNmz+xF4ch
ZKT/nNPtMkzQWGKbiTCE9k2u8dlRDtdvFjzMaw6H0/YMOn+kAvPmNTSiNAJ/
T939MvjwNl4KjOLNUk4h4POH8vdplWCZVNSY7Y+E/Aa3NM67Cm4a+R055nzQ
1Ws+K9beAXs0bbh1lg+KZl5m4sgdjCnnDOPlURiUt0VmNVSDeXRkJ5EdjTbq
g/L0+LtYFrYh+5l3DHK3zazXcK+B0FZ90muLAL/93f3l8r9qIaKFi8+px+K7
H6OM9EdroekTPJv+MBYc/tE0XWkd/tkc8vInSRzaLVc4KpPqEWjiebfgoBBn
C4IOTDk04N+BasV/6MTDWva1NJLSiMuTfhuUA/HQSvF1G3jeiPoklqn13gS0
qGrYtFTeQ3H3c+6rqgTkU6OTJVH3oYi71bVvuQj9Y2n68TukKBF0Uk4Fi6Ae
PPBxVLsJtVbmZ8Yfi+CUOWdk3dOEHCNmx4H1iWBsu92u6f4Aa4ItxlQzEiGe
2bheRfoAaa5P2zmKRLxZml0/aiSDsLKya8g5CbUJBWy9MzIsZM4HSMqTcNHu
nN3WcRnO51RzgxYl47rXXgdfFzmo9a1/rvRLxqZfyzhdtXLc970qmpQloy1L
Ub9hFYng1+E9WgYEnLUmvg4bkCjjiP/KNiTwWDRud2U1iY9HTJXr1hDoC3nb
qUUlUcz8oAoqgVHmi8G3piRiA7tbw00JfDP3y6LirSQU1PjVzywIsI4XHtZ3
IzFkZ/m50JlAb39eXheLhF3XevoGFwKerNzXxAES25KbYhtcCfgzLgZ/8iCh
alw92s0iELHyrOiJN4m6kU/RaocJXOiMui7mkthszwnwP05glX3ke9sQEt//
csDxQwCBvIZwixkeCZgOUQVBBEpKQu4HhZMoWvGJzDxBoE507IlDDAln+7FG
2UkCdKWv3pyAhL7hUs7ecAKyEB+v6jgSkfoDy/siCLR7eQ4bi0hUG8S4TPAJ
uHR7mD1LJFFxas9rQTSBbqY7LyOZBFvSelJDQOBgE6uamUJitYliOjOWQP/W
fdPzp0mMtrQL1ggJ/BdbZufs
       "]]},
     Annotation[#, "Charting`Private`Tag$50800#2"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.885678891644785*^9, 3.8856790875241003`*^9, 3.8857235491935053`*^9, 
   3.8857236410441713`*^9, 3.88572373068569*^9, 3.88572393623405*^9, 
   3.8857515101212463`*^9, 3.885752610571755*^9, 3.88576356328821*^9, {
   3.88576676857614*^9, 3.885766791353014*^9}, 3.8857670586932297`*^9, 
   3.8862454685062437`*^9, 3.88624558599448*^9, 3.8862456568022203`*^9, 
   3.886245724800027*^9, 3.886245809995921*^9, 3.886266173400058*^9, 
   3.88626620585277*^9, 3.886266831516952*^9, 3.886365643114294*^9, 
   3.8863660380591803`*^9, 3.8863660734844913`*^9, 3.886366129855394*^9, 
   3.886366228963622*^9, {3.886366264684868*^9, 3.886366347677162*^9}, 
   3.886366593243637*^9, 3.886368584666953*^9, 3.8863686926764793`*^9, {
   3.8863688168014917`*^9, 3.8863688309465446`*^9}, 3.886368946561965*^9, 
   3.8863689807463303`*^9, {3.886369661642117*^9, 3.886369686271101*^9}, 
   3.886372266202484*^9},
 CellLabel->
  "Out[538]=",ExpressionUUID->"9f030a14-c0a0-4068-940f-74c6a781cb23"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"L", "[", 
      RowBox[{"9.81", ",", "5", ",", "\[Theta]"}], "]"}], ",", 
     RowBox[{"L", "[", 
      RowBox[{"3.72", ",", "5", ",", "\[Theta]"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "1", ",", "10000"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.885678884904017*^9, 3.885678886920765*^9}, {
  3.885723625611165*^9, 3.885723635869562*^9}},
 CellLabel->
  "In[539]:=",ExpressionUUID->"d92a3426-e14d-44e7-be99-e1fa40735cff"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVjH081AcAh9U0fGxYK52MuxMTESZkHff1rpJyeVmU67qYk5cWebnifpeX
nTBpvaySl5pzyijEsHD3q5iGCEtksSWRt9Mkwtofz+f55/k89EORrKCVCgoK
hz7wv3/Jdn7XrTDF7CkO7+6/18Fcci3dynfWwJSJXrfh0hjzDE91pl70BUL0
TJeSjWeZn7lW9giD9SFsbdK0qFhk+v85ZPmzyBh/PmwqaOKtBKtS5M71MUf2
Wy8Xta5VmH1PjY0OtsTFRwGxrTQVMH/w7LogsgVHSDgoxajjn54XKyWXGPBS
VttwxfEziKgJFjU3mBh4l7EpT7YavCc+TeHDjpjIX4iyKluLcl6y/hnl7RBl
c2+yd2hDUe1RHSV7BxSGRbSL4i/gW669t0DLA7prxNGblnWwMF8uLDfejdqe
lFyjvTToPTgaGSpm4cqJxbq+ET0c8cv7I6P5Gxw3UDhd32+Ic0YulI6MfYhP
JQvHYzaifmGUu9bLH++lslkXDSNo5FsvXO0NgIpWs2sWjFE50mpUNhqIN882
S7b+ZIJlRiIZknoIlmrWqYxz5nBvK6GlS0Jgs8Z0j+nMFjxeTQ+dWM0Dx3tr
gaKdFQL9zld4JfCwUTlb43qqFaKfJ7hSWKH4m0YtsNWyRr7cM6xw4Qi8uBnP
X22zwZzmdHWDZyQ8+Cb+1DhbdJVdNuTPH4NhaCN76JEdUuYG6kvz46HTzZyX
qDth6AXlfP2qJCS8/eNYr88OXL5W7TannYY58fCqtNQ9qNEM2+JclQXDX9Iu
yxW88U+3jD04eRYVZr1nquh+oJ70lxW7XsBmPfPMTUf8MfwwjNuXeQn7PyyN
sg7AoOT2XnFODsbeeHgYWh/EupAVHP54Hvwot/Vv1XHg8FQ6MHOqABZRGzvr
aFzY+cCErXQdVEqj5tb4w1ht76zcXvIzkl8mm3BfBUEvcjBd1VqMTb6mXXTG
t/hETTHX43kRHAOnBeniEPy24vBu2aliJMlnPhpaFQqJ3QGJospNdJ9UpA1s
P4KIl31ZO2+VYL6l3UxUGIajvTO67G2lcPRpKpcOh2OnbTYj9mUZrAovU08w
I/F29s77wfu30Gj28SG221EE/ztZOpF2GxsOtiSGDh2F7EaYJCGwHPIHLTfp
ad+BoVj1r1y/Apc4un9lbjgGinLu2VMvK2Ac56T5pOYYBmXNMedqK5F157ll
VUAUmumNpRmCOxCrBQVXTUchx2beVMmnCikq+yjmP0bj8dvO91e+rMZ4ww8j
YvPj2Lg/jrputBrB0nvpObLjCI09mLam/lfAzT5zmh2DFuu1bnNJNWD4fdQt
n45BZn6I96xLLfJZwmseRCzspEuSGEodfJn61q+14qCayvF61l+HuwcO0BMl
cXiwQsnhwe3f4Fv31df2jvHIo8cnn4+7i6TDruVrOuLRN5a2TrCtHg95Mdei
9/OhGP5selStAb+XmseOyvnYeXaRatfVgIYbn5SXJZ4A0+ZWi7JPI77ULbU+
vvYk0ufNTRXqG5FkEKzSKjmJYfWfakapUoQP9JhJLBJQTeTzNE9L4VlGFYTd
T8BFpywnq3EpOqZq3do8E1EcsMuF4ynDd2H9nbNdibBoLwltq5bh1f6+BnmQ
AM3n5DWb15Owd20YrnotgIfq1NILbRIqe+NE6pMCPBKOO13VIRH/OZcWMi1A
b8RIqyqdxK4nHV9pzQow6j4wOGJIQj/FYoi/LMCni7+rXLf60JdOWemqEGAF
Fexb50Vi/bK2Vr82gZ6+3Nw2FgmTe5+3KekQ8Gfl/J3iTeJp95ZES10CXObF
8Dd+JCqG37SfphGI1soUdgSSKJjg29oaELjQGlecHkaC16RucH4zgfXOMROO
ESTKlHszGs0I5NZGWc5Hkpgqfjw5Zk6gqCjibkgUiYFP2UWOlgR+FR7ucOGT
0DjY3z9pQ4Axx9FcPEEiT2yvq21LQBrBDqhMIHFmV4+/69cEWgL8X+gJSbyu
LbqfwyDg2eln/PQUiR3E0EyTHYFOd5/I7GQSY9e9dWbsCfg2sCrdU0moFi45
6IJAn9Wed8vfk2hfOcjZ7kDgPz6g4kA=
       "]]},
     Annotation[#, "Charting`Private`Tag$50849#1"]& ], 
    TagBox[
     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVz3s4lIkCx3FJKCesg9BmGkehyz5KZtuD5jcG6bJkyiWsWZSMYVTkVtIM
xvBSOAwbobOKOhaNqUnsYN5CtaUQR8qTPbk0IbHkWmf3j+/z+ftLDYpkHVdV
UVEJ+rO//CXHaf6FyiR99vSqqBcWA/TPLtW7E5x08bJasqGj+RM9m6M1LRd9
DVud/BTDDlV85SLt4YeYo5cpj6qwWQ3f3t9tykVb0KZ6ouAOXQcsqcg12NMa
jNblnXdT9TC7RImNDrFBfZB1rpBlAPpFt26x6DtQTxSX2kpM8LZnSLXyJ3tQ
R4j7rfs3QERJ3FF/kw6rialT53oo4PzXsy1i2BGV6jvVo0PMIOGkmGdr7kMR
WzeVYWkBNe1nDUY5+9Hd/6/fGtwt4SVZf/iq8UFMCyjzlGgrLC5I+JIt7iDm
+N9zpVth1noyMuw6C3E8ObNWyxpc79LfMtt9wLW8uSmTvgt5Vs5GzzOP4hD/
5vZAxS7IF5XBBh6+0NccmkpysoVuGW3xSp8f3mQ98q9zoUE6+sSqRhkA/wNb
ozKYu/HF/jwZKgzCfsbiwksze7g+rdpIVIbC3K/DY7iVgS49atiEHgc1i+3l
F1c5IsA7v84jkQND76Q6dSdHRL9JdDFiheEHZreE0+yIsim38GuLXJwd2+X8
dxkTc4YfZU1ukQgwcXhVUuyM7prLFgkLp7HbY0Vi38F9SJ0bkFeXxYPd8VDv
yrZD+H3IKF++KhlDgc/17/t44fK/ZXvn1qdj7lQy77KWP+oNw3c53bmEJbo7
x2/+R7x9oWAPfshFnntzR29ZMCjnfBU3XMSIPCpN6w8LwfDj8OD+rJ9Q2Dj4
VN+Ag01Vtw5fLy6GR++MpXUyF+tCVwQmjJdCPLRWwDTigfGyZWBacBWndJ5N
GNudhIMntrE1fsaJw/W1HPYp6O1x0uyoKseaNRvVOwtOwyxykNCiXUeY/pkm
R1kU/qatVnLwTQVircSd7cpoNK445q4Q3ADP+/I/3unHoNLhh0q11f9B157I
wi6LWPBG+i8dqK1CEq3v7Qu3OJzsmzZl21XDx3Ttg1B2PA58l2MfO1IDl7XL
3LHkBHyavb00+KAWurUKfU7ZWYTMfKieSL8FeF0ViCrOQXEzvDIxQIJb1Lg/
YhSJsFe7MzNlXgfmmJ1E1H4eRpoluYKROsy4tk2K25IwqGiPybsnxYSmvHBT
xAW0U5urM5NuY3rA4V3E2AUUf7uwXcPzDnyE8V2nz/DR9alzqWizDDaN4iDd
RT4s/eMo65Qy7Mk07lbGCBAW+2O6vvwuTNX21nGXBHhEM9g7l1yPTUXjqw7F
JSOrLPTIrPM9DBTMZZp8TIZDy+fKGKMGODePPjwTlQItYaDH61cN8P9lfCJy
LAWtKzQYrbcaIX687f4wLxWl1PiU/LhfcWV0h/760VT0v09fl2Qnxz4ajcoO
EEIt4vVHpXYTXlksO2oNCHEgd5ni0N2EJwx/5rUjaaB/W/tI07MZNLOglfnd
aSAWrLeryJux0e7xlyIPEYZ1CuuVlBY0SFe+b3wgguxCGccwowW3X9+eoTDT
UcC8xLQdb4E6fw2lTpaOG37fOwe6KUD9eFwjwyoDOzqqwp7KFNh+rigvuzwD
7XlT9d+YkJi35aa9NyRwUGvy89B6Ev/k2u94Z0TgGX+ceWUDCTdhQ+ewCYE+
3ugTLSqJ0pnWpUFTAkrXgcFRCxKTx7aO924msHb54eqfbUm0e32gK2gEWMev
Hl3nQUKgx6/L9iLQ019S8pRFojux5UiWDwFfVvH/Uo+Q0Kp/PJ7uSyCYXhDx
hzeJr79iaCQHEIg2zuI/DyCh7d/WGx1CQPwk7gYRTkInWnXSK4aAiVPMhCOP
xMo3LtqH4wiU3IuyWYgkcfxZkLl7AoGKCt6voVEkeq6Z7957nsBd/rHnzgkk
5OMjyzQhAfu5QMPlsyRqXFX6d4oItPDYftJEEplTszXfZBB45Oc7ZMYnYZov
YGy+SMCt03vLS8GfP3uMP1OzCXS6ekbmpJAwK8yRbMgl4NXEkroKSSjjxwKM
8wj02x6a/5JGorxhy0oDMYH/A2fG6EQ=
       "]]},
     Annotation[#, "Charting`Private`Tag$50849#2"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.8856788916866617`*^9, 3.885679006140519*^9, 3.885679087557064*^9, 
   3.885723549225946*^9, 3.8857236410924387`*^9, 3.885723730729806*^9, 
   3.885723936276074*^9, 3.88575151015632*^9, 3.885752610604838*^9, 
   3.885763563322987*^9, {3.885766768620377*^9, 3.8857667913925257`*^9}, 
   3.885767058739169*^9, 3.886245468547699*^9, 3.8862455860315943`*^9, 
   3.886245656837278*^9, 3.886245724837285*^9, 3.886245810038374*^9, 
   3.886266173457501*^9, 3.886266205890098*^9, 3.8862668315659533`*^9, 
   3.886365643156186*^9, 3.886366038093513*^9, 3.886366073522078*^9, 
   3.88636612989297*^9, 3.886366228997745*^9, {3.886366264720487*^9, 
   3.8863663477123528`*^9}, 3.886366593282444*^9, 3.886368584706367*^9, 
   3.886368692710593*^9, {3.886368816836722*^9, 3.8863688309812527`*^9}, 
   3.8863689466057777`*^9, 3.8863689807892942`*^9, {3.8863696616824923`*^9, 
   3.886369686311803*^9}, 3.886372266245026*^9},
 CellLabel->
  "Out[539]=",ExpressionUUID->"b952f162-66be-4bc4-b8c1-82bab1839b72"]
}, Open  ]],

Cell["\<\
Trying to work out what the implied viscosity is from the assumption of 0.38 \
by Wilson and Head.\
\>", "Text",
 CellChangeTimes->{{3.885758619286591*^9, 
  3.885758648934758*^9}},ExpressionUUID->"41a96361-fc1b-453f-89b4-\
e498b0e27f52"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Tau]", "=", "2000"}], ";", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"2000", "/", "\[Eta]"}], ")"}], 
     RowBox[{"2", "/", "11"}]], "=", "0.38"}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.88576596661662*^9, 3.88576597028648*^9}},
 CellLabel->
  "In[540]:=",ExpressionUUID->"26b42d93-c524-4eaa-8b86-9118817303be"],

Cell[BoxData[
 TemplateBox[{
  "Set", "write", 
   "\"Tag \\!\\(\\*RowBox[{\\\"Power\\\"}]\\) in \
\\!\\(\\*SuperscriptBox[RowBox[{\\\"(\\\", FractionBox[\\\"2000\\\", \\\"\
\[Eta]\\\"], \\\")\\\"}], RowBox[{\\\"2\\\", \\\"/\\\", \\\"11\\\"}]]\\) is \
Protected.\"", 2, 540, 42, 20715073998163303067, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.885763563330962*^9, 3.885765970822875*^9, {3.8857667686297197`*^9, 
   3.885766791402169*^9}, 3.885767058747472*^9, 3.88624546856988*^9, 
   3.886245586040084*^9, 3.886245656844408*^9, 3.886245724844943*^9, 
   3.8862458100460663`*^9, 3.8862661734881783`*^9, 3.886266205898209*^9, 
   3.886266831576441*^9, 3.8863656431758757`*^9, 3.886366038102407*^9, 
   3.886366073529155*^9, 3.8863661299011307`*^9, 3.8863662290046873`*^9, {
   3.88636626472759*^9, 3.886366347718893*^9}, 3.8863665932890587`*^9, 
   3.8863685847141857`*^9, 3.886368692719296*^9, {3.886368816843638*^9, 
   3.886368830987685*^9}, 3.886368946615244*^9, 3.886368980810235*^9, {
   3.886369661690873*^9, 3.886369686322673*^9}, 3.886372266257203*^9},
 CellLabel->
  "During evaluation of \
In[540]:=",ExpressionUUID->"fa08f2ba-3efe-4e3c-8b8d-182fb14a497c"]
}, Open  ]],

Cell["\<\
I worked out with help from WolframAlpha that nu = 409468 is the approx \
answer
This is ridiculously big. Why did they choose that? I think it is to fit \
Figure 7 of Pinkerton which says this stays fairly constant for terrestrial \
lavas.

OK, I will now try to include the nu viscosity - note this is plastic \
viscosity\
\>", "Text",
 CellChangeTimes->{{3.885758553806971*^9, 3.8857586138137493`*^9}, {
  3.885758651846902*^9, 3.885758658574898*^9}, {3.885763379942031*^9, 
  3.885763418639242*^9}, {3.885766356925324*^9, 3.88576639544461*^9}, {
  3.8857668395165873`*^9, 
  3.885766846867885*^9}},ExpressionUUID->"4870c497-cf11-49b7-9695-\
f62d4fa430c4"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.885723671635119*^9, 3.885723671636828*^9}, {
  3.885766231161447*^9, 3.885766247321267*^9}},
 CellLabel->
  "In[541]:=",ExpressionUUID->"6fbb851b-6568-4c96-9532-ba8ea625b693"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Lwithnu", "[", 
   RowBox[{"g_", ",", "\[Alpha]_", ",", " ", "E_", ",", "\[Eta]_"}], " ", 
   "]"}], " ", ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "/", "1000"}], ")"}], 
   RowBox[{"(", 
    RowBox[{"1.34", "/", 
     RowBox[{"(", 
      RowBox[{"0.00007", "*", "3"}], ")"}]}], ")"}], "*", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"2000", "/", "\[Eta]"}], ")"}], 
    RowBox[{"2", "/", "11"}]], "*", 
   SuperscriptBox["E", 
    RowBox[{"9", "/", "11"}]], "*", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"2000", "/", 
      RowBox[{"(", 
       RowBox[{"3200", "*", "g"}], ")"}]}], ")"}], 
    RowBox[{"6", "/", "11"}]], " ", "*", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Sin", "[", " ", 
      RowBox[{"\[Alpha]", " ", "Degree"}], "]"}], ")"}], 
    RowBox[{"3", "/", "11"}]]}]}]], "Input",
 CellChangeTimes->{{3.885679102555542*^9, 3.8856791064013433`*^9}, 
   3.885723723777371*^9, {3.885766254082788*^9, 3.885766301612267*^9}, 
   3.885766347070053*^9},
 CellLabel->
  "In[542]:=",ExpressionUUID->"bccf9706-289c-48df-8278-b45694d34d25"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Lwithnu", "[", 
      RowBox[{"9.81", ",", "5", ",", "\[Theta]", ",", "270.43"}], "]"}], ",", 
     
     RowBox[{"Lwithnu", "[", 
      RowBox[{"3.72", ",", "5", ",", "\[Theta]", ",", "270.43"}], "]"}], ",", 
     " ", 
     RowBox[{"Lwithnu", "[", 
      RowBox[{"9.81", ",", "5", ",", "\[Theta]", ",", "492.75"}], "]"}], ",", 
     
     RowBox[{"Lwithnu", "[", 
      RowBox[{"3.72", ",", "5", ",", "\[Theta]", ",", "492.75"}], "]"}], ",", 
     " ", 
     RowBox[{"Lwithnu", "[", 
      RowBox[{"9.81", ",", "5", ",", "\[Theta]", ",", "665.14"}], "]"}], ",", 
     
     RowBox[{"Lwithnu", "[", 
      RowBox[{"3.72", ",", "5", ",", "\[Theta]", ",", "665.14"}], "]"}]}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "1", ",", "10000"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.885766736317161*^9, 3.8857667866673594`*^9}, {
  3.885766902575794*^9, 3.885766906494794*^9}, {3.8857670201215887`*^9, 
  3.885767054434577*^9}},
 CellLabel->
  "In[543]:=",ExpressionUUID->"3dfd632f-5932-4054-979a-5d6182690d97"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVkXk81AkAxfGxzmpTOdZkGdW2W6uSWGrwbI4pIsOkpVWuFmHs5m5VzkYo
rDt+rsq5DMYxFcb8PpjaJQq70hZyJeTY0Cit/eN93l/v+/m896iuLIaHhJiY
mOu6/vffk03f94nNGQe+jBd19cwZr5lX6oeZboay7PJBkytySPKSX2xmb8dt
V0Wh7eo2KJhz+yPO7wRfQ5w1t307HP8a0bnN3oNv+c8XXvtSweCy6W7MA+C+
0hKElu7E0gf14IDzOnjPpMkpPdwN4xvWvelsA/xcLfm4qkMLo/1jEiVZNGTN
S4lHnToAtnq4Nq/MGB6MdyYSAm14/c3s8B3/Ht2N1H8zYw+hxit6Z5LMMcRZ
7tDe1mcAyU3d91WSj6PQ4O8EzpdHcKqGYlfwhRV4kGQPutOwKqqJqNljgxTd
tIu8CSNotvuzvO8y0NYr0akzboILDnl/JghPo395VOCQbo7Ub8xUehJ+gLuO
3OCxry3QvDrlpmjriKKioEJ/ngU25+ut5g444eW0FD94gA7uZOc3VVPOSAgr
yldSsMQn2mXSM9YVsoo3pN+42oDeVaERX+IJiph0T06vHZ5uoXrPbvFC1nUF
lUhVezg7pNXahntBITTkmvY5ewQMhZurMLxRlbkoWpiyR/6Ctc+d1Qu4pDQg
TVljYkVpvqHFmoV2WZPoJspp9FZl7w4T/YKkuvPJ8gZnELPyorkyPxQdfiPH
Z+VcMTKmktb8WRQ2FKbVETOeyC5ssFihxKF6SeKZiM4CT8nnkGn9TbyZVQuw
e/sLRvsEZ4ffpmBcGKBqPRYI9V8dBaXm6cjYIaXVKR6K8T983AYTs2DpvFzW
J7qEXRXVdndzchBoI/Tq4lyGsqe4S9hMHhI9LgvVH16FybPWF4uRBYjpXM1Q
pkTCkIlvz0oX4WNBy597XaKwxchU5nHFbWTrU4fKmqKhyRqOl9e7i3j2jjWb
tRhs2CRJWA0Vg6LLcTnDvIYH4u42gshSaNBPqM4UsFFi+GOJpGw5tqu3JGf1
xsFvYvCmJacCuu7aShy1ePgPLH559kglCg921LWcSIClQTIteKIKaiq0ro2p
iVheqvsw3MZBMfdR9E89N3D+3dvK2bhqzBlcODS9dBOCMp+ScOcaaPgODZhm
JYEmWf9uYWct+FZv2/z3JkNFhkiJnKiF4lyUYS0vGcMCYVDqPS5y3GoN3zBS
IKTyKxOu1OGwotqdiaEU5Hwn0pJm1oP63CiDH/obni4/+XDrqwZ8khyN6dyU
iq/PhKgrTzVAYavFnaXMVHgHn4vb1twIj+P5v73SSsMjPUWLlSgejlacbFRu
TENivqf9ktk96HRznfuN0mHYulYSpHIfthtkstGdDvlYF9t/nt9HyHSS8yGH
DLSLS5u0Vz+A3DsT16bJDORRQ6PTQppQHiSSXWBlYvBNnPKVI82IyTRlv1zM
hKTvP/NTm1qQ59lHt4xc/z3lo7phbwu2BpqVZEplw/g7ziMZJh/iJx/026Rl
I150QEusmY9o1oi+UOMWxj/P5E2pt2J/vbaLXeEtNFzN91K63go9zZWLB/fl
IOPozaO6M60o5zzbYMjJQanTCTMXawGsshyvHt6fC+3HFd5dDQLEluibBDfm
Qpi6wNunSoIMW/0krU/ASn5ubYxCIqJ9Ng6HCXRHzBzNVSOxayJLLpRGYMBv
slOeSkIvgzL3GgSm6C+GJ3eT0NSftPzjGIGNHx/KFumSmJaNtEhwIsDwKPhB
2ZbEE9vrTOkrBPoHCaKLQeL1+JlGowgCjoycVzH2JIIDpLYERRFwM87w/deB
RP/JkdLRawQCvkiM6HEmMT/114XWJALpnSGl8T4kKNNJBqGFBFRNg2a/9yPh
MG8mx7lNgLh3UUfEIjGiP9wzfpdAcbFfk+fF9bxjl5l9OYHGCPceszASTwNy
+fu4BGgrLkofL5FoU86je9QTaPU768QNX+cl3ei41UjgkZPjmGbE+h7Cw2Uy
TQSsnzjseRZJgu/9fqNxC4EndCYrOZrEeHGpZ2ArgVMtDC49dr2vv9W9cpLA
oO7J95+ukdCvH5EYaSPwH8Wc9vU=
       "]]},
     Annotation[#, "Charting`Private`Tag$50906#1"]& ], 
    TagBox[
     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVzXk81AkDx3HJbjS7xMZMVLPjYZUnbbYki3yLkWJZoyGmpaHELKNFaqdV
uUpR7iuMcY35TY5CGMeI31N5PKntcORoqdybLG25d/f54/N6//lheAexTigq
KCh4/9P/LUuyme9UmLby9W52cI9bslqxLd8jsFmHtBvVISPntJDoT5mVx27E
vZU5zccdDKjbVndF+OphvGuz4VSvATy6X+0sijVE6nQmzTbVCKzqWDsf9g4U
qDu5HY01xocl+plQ351w4jd0D2zfBavrjs/TY82gPOQUS2Ga4U3XsKIkywJ+
AQNHtevNEUsPN5ZJrbBKzSdzWXcv/HvYDwJH9uPb5yo9JG8fKv2j9RKVD8J+
QZZwYZIJJdVfG2hJh+B/wkt3x+oDcK3Uccnf4IC1NvFP3lPtsLhQGVFp6ITV
y+9sT5sdgu79U0E8MQtG3tJ1vGOO+NEt72F82xHo31OUpRq6IHUrk/Yk3h29
v3VuCRC7QL444aPp7IGp3FmNDMZhrBPtXsx9wYH73qTSFU02qsc6tlZMeKKt
6WJi8UdX/GVxnvS75A3BRPuRlXwP2D0q/TJO4gcdUUtW47NjeKbB4E1p+KNr
LaufosaFp1talXO4P6Lmu/hGh7gIHQy3pbF4cKlUV5hs5kI04xhQvPgjpncF
p/sS3pjT+qO22TEIrgqf5SmcPo7nFTcMBAvBSGg3S3wwcxIxcy/l5aKfYTw1
OyMX8vFqmJYm/yQKZXLdgzSD07hRUHtgTucKOsTJvYtxAsi0AnbZ1CQgy6TJ
1Ej9It50tnoNvUtGoemYOSc1EvRfPFoJ23Q8NYi02bMnBiP/C/Dpu5YF7UFe
zXLTZeiX3nYR5+RgbFw3rm/tVVD9VnEFb/Ng3VjBzvklHvt6W17ORuaj4uoR
l23112HJxjavNYX4Ou5czIOfEqGx10b5cWkRxjoDk89nJUE3aCiOsluM0sfM
yylEMj5TVRI6DJbAtKKMfvNhChpXHXdqjSRAG9gfMTCSConlDxIllZvQmKyK
pE6mgT/al2B/qxT94/5RrooZOPVidrOXeTmix4TzR5QyYW+WZHFmtALGb3NC
zTdk4eOHO0tD927BwOE/8uivbsD3z3flU1duQ2XUUCWNno1WaYAk3LMSnVt9
fE78OwcWSjV/zuhVYVNGnm/Yv3JBUxYmR45WwY16qvf1eiGGWtvCUuurwWZv
vzZBy0Mb4255/IU7+G20O35aSYQc0wWjNewacF8SloWXRHj28elS9le1mBBC
d9Mn+dhy9CydOlGLh9kxP7kL8sE7c+zKenkdrqukbK+by0f7bs0Dc1EyqHX2
dC8HFuCayO/wB2Y9xDnn6IcHC2DZsiIJozWgcO4b2zpOISiXuM4D/Q1wDROI
zB8V4v6qNfvu327El0zKiOxgEfIYP0ennW3CGuFNp9DGIvRNXqFeMJdDoq0v
8d1WDKXAgT8mVJvhquPMQ3Ex7JOX6ZbPm6HOObM5WlMMK9Nb7crsuxDt+DiU
e12MuIUdRgryuzC9o1/2QakEI2qZsgl6C0xKH9JfB5eg9qLIX+tqC+LTv1BN
GS9BhnWCtcnbFowyDYsm3SQgON8xuY6tiFVt8ikjJTB+XMp7VNuKnYfWMxh7
CLSlzsi2a5MYki30exURcKBMrwzrkPhhtZv9fjGBXyPeWuduIsFd/ey2noTA
C/5YB4VBgn/e6OT4TQITdi+HxgxIfN0QFRJcReDz5f+qFJqQyMwO+iKqlQDr
RL471ZnEi37VMOEgga4+ofARi4ShxfL4hVcEPFg5r2MOk2DK9V25bwj4WGUE
vncjURy4hao3RiB0w7WIJ54klqwpFsQ7AukdZ4m4ABInC5bUK1cIaNuETe3n
k7hcMaWSoiCFsD5k50IQiageg4UQRSlKSvhNfiEkemqCSZNPpaiLOP6EKSDx
+9aI+brPpbCY42otnyMxnt1alaUmRQvfi1MdToKT4ewrUJeineMxrBvxz8+T
U22uKYXjUzfD3kgSiWS300aqFE/t2EFJ0SSshoteL9GkcG1mVdtdItGbQAYN
aEvRZ/L9/F+XSRQk7HrftFGKvwH50PTo
       "]]},
     Annotation[#, "Charting`Private`Tag$50906#2"]& ], 
    TagBox[
     {RGBColor[0.560181, 0.691569, 0.194885], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVjn081IcDx1kWm4dJ5blElti0mjgN84k7tKScHTl5fhi60JjHqY5LHipO
OO7r4WITMYRCcrr7/ob8ykMeVmk2lYcpzxGi1v54v96v939vbZ9Quv9HEhIS
Ph/4z79xqSsDErOWK1Ktz/P1py3f2VSZxlIVUXKC8tRYLI2MINkFYbImOJvH
c77uUcImm/pBdoAu/Jg7N+saq4H5xzOjX5INcI+2K+3yQS3Q65PtfBl7UW35
zVx3kg6W1rSiIgKMcGyZye9w+hyWlx36c5IPwNBc2NNRZ4AXg6MfleWZ4/5Q
iQHb3hDJWvH7mq5bYmCh9zv3R18h6BGj/dSYFdQ8KK9SfvgatUEc3QyZQ4Bx
pnWRPgVSCj3NqtzvENFWvWLgaArnWg2nq2r24KaaBgh+OoC3q7XsWoOjWNXw
uS1/yww6bWGhwaV02Bnu5xfLASddiu5f7DiOYKWB1y2gIkufptp70RWa7zvr
+SQVwreTvlsdmWgL+kZPRKNBUWDytuCxGygb79zfYGeD+okH+tWTHmCc1p5Y
ptrhvfkZMjDJB7M9bS8Tde1h11W5I60sEL5FKl3z7Y7oU9IOnlYKQv9Jn77C
jXR4uGTXOcYHwXM1g6tOoyPi73gbVXowjCQSxzkiOgTzDqxf357ExCMzaYtG
JywrzzW0OoRieO4n69ECBvqr+Xqxqz+C4T6V6+LgivPLw8IqQQzqZhKZXEMv
PBtVzRZ+nAjRRn72VmYA+MUNtssaKVjga0tYybPQpMzaT72Vju3Ftr+7rYbh
xYDYc2QmE7zNklYSxeHQ+pkpLrfJwd7uJrU9rEiM/Z/lO3QpD61zyYJ05Rh8
XnnDqTQ/H2E8+6puThxUAiW9Y6eKMCQqkv5U/QwOPhENLyRcRd4O102fmZ+D
BQNfekqXoERPVinMiw2lb6ky3ZW/oMGkZWAyNwE6oSNpsialIBQOcViNiZBT
kCq0//satvTMN0u94uCOpN9RcUI5WmK4JSbKSSizcC+T+qQCtBcRx/bpX0DI
+FD64ZpKxA77SdGOJSPs8cJ2T7MqdPbYrox7peDwAa551Hg1frNRfPbz+VS8
Wbq5NvJ7Ddqfnnb6tDgNAYszVdMpNxAmsLI9UH4R4uussniPWmzSoQSb/u8S
zKVuLc7r1qF5dkOha+dlqMoUZiaM12HzP9cnIu+lY0TcEZl1ux7THduqt4Vk
oEP7btXFszcxMbjEPjWVgXzKqqE04xa84/LYsZFc9L15uEbsakBo9s15rTUu
dp+I1lKZbEDD9vZayehMBEd5pWwRNuIgxerH1PVMdJpstV1ObMKuFrktkTFX
cEkQ+P0S7TYGFwOabeevwEL0rixStRnHsa38ekQWZJO8Hf982oyWTMxVTGWh
TVL6YNuNO/DljUUbhGWjSDuGkx3dgtL86Bn/f7Ix9DJF5ayZEOs/lC3VeOZA
6tSfc5MKrZjSrBT4/ZWDw5nrWhb9rTjZayf42JkHS0pNpwzjLiRrk1+tDfCQ
trrXUEJ4FxoFZodknXIx9llu06SWCEcsjisYteei4ZwgSDn1Q6eKzxRQ88Cz
Trc2nhIh/k5jMqUpD+VuR2jeDmL0d/rv3voFH/u6K4O7GsQYfVL+RqeUj46s
+aY96iSagrLLOaoE7GVn341qkDiqFH06QZ1AD3vKumAbCfmdMvvOaRJ4HDLx
QFabxJjUfl7cDgKTdsMjE3ok0qP63p/eTUB+/d4nJcYkul8vWrmbEqD7X3VV
cSSRYjKeY3icwOBQYWEXnUR4jL6zAZMAk57//Pz3JCwVjeX1ThDwteSdeu1C
gl9R4bPDi0CE2iV2rwcJzeo9FUqBBHIeRJensUiclbWvW4wioE6NnLYKIdGZ
lTozF0Og8Ha40WooiSN6gp3TcQSuXQtpCQwncbnbKXzsLIFGtl8vLZZE6/Lr
msELBMyXvZXX40g0526ofJhCQBTi6VYfT6JPbrawK41ApxtzVIdNQtDnz2pL
J+Dw0MXgScKHnz9mbcRcAg/tGKFcDgkTRXc14RUCzq30erskEi+dK543ZRMY
Mj628v4CCbnEp7/e5BH4F7Xd7rc=
       "]]},
     Annotation[#, "Charting`Private`Tag$50906#3"]& ], 
    TagBox[
     {RGBColor[0.922526, 0.385626, 0.209179], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVi3k81AkDhymVaFVqmbSt4+1Y3lUkVMz65ihWKVNjlinjKGsQfdQS7ZQj
GoumkYyc45rfVI5I7hjz22S960jSMa3SRrOK3I2rtveP5/N8nj8efZ8Q2olF
SkpKPl/4v0v49jOPlUZtJtNa/vhP+ozNp72lOyPtVyGQuhSqhmtxha0+0cj9
BuX8fUnxw7pYvbeyN9pvI56ORWW3+m6Cx5PXZoVcI9Q46DHL4oxAq+Q6+tJN
0NqhGJg7tQ3T87rhZ/zMsMe7bOwlZTtsLrv0pHF3IeyW5FiE2AJvegcWia9b
Y/y95eTFFzvB1eWY1t60wXKp86hoiRXYT+kPTg7awirAd1JuZYMK9sWNV1Sd
0JZh+CKZsIOKRlc9hf8j1jdY8t3y7eFWsf5w3rr9iK36VBKV4YC52YroCqOD
yBJTrgpj98Gg5VRIgIiGrsUy5zprZwQycv9Mav0JtHNCF842V6QaOlAeJrnD
7PcSS26OKxrnhny/dvXAXF9siGIFDauEFnPZz5j4QZJWs0FOQ6W83bBsyBPb
eZnuC4Ij+Gx9nvSP98Hyz35yZj8Djh3Feolif6i1QkfT2BOPNPUDRjTZWLCY
tgsWeMKTce2OK4eNmCaDC2eVWTjzirOXQgvA1AWtgOweFoTjLkFFc4EwS9cb
F9d7QaE1Vt3kEoJyrvY3JMMHPWUZWyJnQyHdnLr0vxEnEKfoaywVRiDAx31F
zplAvB6gXGtcEouvmIPy0I+hyMiv3qdYn4BZFfuG3PZw1GoF7bCv4kFtlSTU
qOtXvHksZfV/SEHmMmVWmzgKur96SG/sTUPM4ojO9qWxGPxfkK8s+ToW21G7
HzvFYVNx+WFRVhbW/Kl3/4n7JWj7K3tHDufiQKGg+f1vCdjzvLlvIiYPYdOH
7psnJ4JKx/esZQVg35K4/cxJhuYP9qqdxYU4apdZf5fHg0FIf6K6hQiKA6PU
h/VXsEJDJWf/KwJR9+Z+fCTho0H5+EFpzA14SXwjZfdSIKYeE6ssvwWd7zSu
mBBXEfxWxnO+XYyZ3rGQq3mpOPVs4luWVSlWV27If8e9BuddfOvwt2Vo4XmH
mUSn4eP03fn++7dRLI41KAoVwG/qQ+lIQjkObX359KNrOqQ3g8QczwpIzz9P
0LG9DmuVqqnxjXfw3kmNydXLAEU1JyXm7R0IakjbRTMZ6Je2hqXWVWLFCVOG
4etMtOpLSpMu3IVsCSdvZ10WsixnjZfRq/Dudu+qzvxsPPrYPZ+5uRrpr5p8
c+Jy8N3Rs7raQ9XYTjsZI2PkIiDcK2FtYw1+4vNXUzcI0Wbx9T5FbC1os5zF
x+OFSBb6H5l2qEP3Wh+rf94LQW3+JA6j1OPl7Ln7Scw8qMd7u/71oh6DpjE/
75bmoUV52Z6W8gbY9HodmzLJR65+xMVrZ+9B+bKyZ48gH7J3CdoXrBpRNDL8
Zsd8PlRO/jU2pNEEt6Omtnz/AjinLOhSe5rwzMy+YXd7AWwsb7ep0iV4YhjK
6dpViMRZE2OlRgkGM+b1RwoKMbgyvXZItxm8Plen80uLUB0lZGv91ozcB9yV
Q6eKILDj2ZkPNyOvlWose1SEG8wDDt4uUmxy8Z+K2CGCaWdxQEe1FIwI75I3
WSK0po7XbtUhYfT7mksTCyLsVx/9NLCexOVi3rrCzyJ0RQ/bZW8gYVEyV3hE
mcCzYHm7uj6Jw9Mi0V0VAkOOff3yLSQC01n0cHUCXy38sbzAnMTpLUreMxQC
tBN57tquJAq2brNQmBHoleXkdNBIaPz9rZwwJ+BBy/o77ggJuZfWVYYlAV8b
wclJxpe/SLenajeBM+uSox96kmhqLlv0yx4Cae1nbyQGkZg8yNYcdSGgYx82
YhtMolvVmJ57iEBO3Wmz2RASU9+r8F1oBAgi+J7/aRJZGR8miukEaqKPP3SI
JMGz7vRjHyVgrfDWWjhHoiraNpTiSaA5mMWs5JC4OdUR9oBFoI3pMWAQTULN
ySJwky8Bl26G0fMYEhUCTWbPcQLdjvQQ/kUSKRNq9rF+BNyaaJWO8SRG1+ht
3u5PQGZ+aObzpS9dRlPuZxP4F41g+Io=
       "]]},
     Annotation[#, "Charting`Private`Tag$50906#4"]& ], 
    TagBox[
     {RGBColor[0.528488, 0.470624, 0.701351], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVjnk41AkDx4dHKGUXxSCJtdvSgVcoOb4x7BQrZs262rE5epkwthjT1OS2
7OTMYGbliDTaiUiOMDK/tmyrHJW2lE3lWBvpIAb1ev/4Pp9/Pp/n+RoFs2hh
iiQSKXhl/+elPMrCA9KM0yfT4p9acl45fXSr3cWlfA7R/Zwl3wRl5EaovZNm
bIR95afNmX0a0HBrHEw6bAJHwcs4wbAOAh4+t6rKMIPlNWbfnWID0BozqCF0
C3yRxu5U4xthbskwPvawFRgMW7NZSxM4ZXveL8zYjR2txfsmqaZ4OTiqKBba
4xGTiPOXbkWGIc+y9aITeijtW3/dsgMRf9FvRY05g/zsXnN0lCUaIlJNclX3
YcZwSf542hpK6n1t5Lz96NDmhTQp2+L7Bv3vKnQ9oPO3PO68/i4syhuSGswO
ILpezj7uYAfjmzEsZjUNukGUtRphjjjiW9ZzutsPihZHdqRud0GBqSu5/7Q/
zBZaBdE1LpAuToZs8A7Ae7fxgxITCj4vt1k8+ygQp2rfhu7WdUXjxB3TukkG
emcN1pktuuGT/SkiPD0Yz/t/5GtW7wf1rmQzXxyOA15f+nkNeuGephFzWjMC
N7ZpXRZreIPhK7jizYtAp8uc1z0Pb8Q+47mRaUxw9TiT52TeKH/rGXl+8Qh6
EtVdnCU0zGu/ae70ZMF46mCdMccH9+tEW7jyo1jv8lnK0TlfpM0PS2vLj+OJ
wn+5rEoGno+SBdJVKVDdNldGNwuD6FzzN/P6meC329auzWWiVTtyJ6UpB3Z/
fDCt1IrByweyoJHX+RjwjaqzKjoKw5MBshq3QjTQu+w4e+Iw9mdkyFCWEBk9
OWfKu+LxpaT+u+qSEihP6Bb/pc6FTrjCIe5UGXaurt+1kHgSex93Db9LrgB1
n4n6iPQUHOjYFqRSic4yexeP2ERoOlJUeyVVCLb1F1mUJMGYNcJXs6nGgP4N
O5YkGWvVlUo9nl3Agf3a4Vd6U9CuEHpAllyD7MwBBY3JVIgdfhArrf4Nq9KU
h49OpyF6fCjH/bIEzu5WlU9W/YyYR+82Be2pRXuY01ll1Qy4786zjx+vw/ms
TYlaBpn4MHd1aeT3yzDIdFBNN/sFh2df105n1oNPao9v/4IP2cVIMY/RAH8b
y8Gr5qdhr9Q0+9bkCrJkkQ9fbckCWbU0P3n8Cqr83lASdLMxIutmF1xrxIbR
wtuXDHLQbXS99nTC1ZX/7mWhKrkosZVvV6E3oSnFZvvWzFzc+zCw9OtXzShI
KIhpUcnD1wc5hjqTzVi+4dOozMsDM/7HzPXSFtTcIvVwFvNw22bDN/MprRhL
sjjSF5OPrPJwnznXa0itqrHe9SIfDl0fxWxyG2Yv7HSuY5yBWvoh76dP2jBY
u/srSv8Z3FRQ2Xuzvh3De0vpDzwKUGZ0PFXA6YDfC6Z3RWcBhv7N1EnYI4UP
T7G3ylwApainbybVO1G37rBWnlgA9/xlQ4f7nVjjuHlhlFwIJ9vLt1Xp1xGj
d56qmV8IvtxiO0l6HZK2Z1Y81SKMfVbcOmnYBfIJn5NcdhGaE8sjtH/pguc/
pGL7qSIUueS4WE91odU8dn11YDFqAr91PeQpQ2763z6cW8Ww7JUw7zbLwHeS
jy/vEaK74G3rDj0CaQVPbb8VC+GhNvNxVJ/AseqNM+YXhehLmnI5a0DAwI5b
rikR4lH0xB01IwJnHAv+eVgnxCR1eGRiC4Gv9QINgpuFWLf8x+pKawL9qdNq
cTeFoIVV+Ot4E3j6em6N4KUQg0OlpXdpBEgdFd3sMSECaCUv0nxW+pNpp/wm
hAhxKop670vgudbigP4rIWJ1s5L6GQQY8Vc3nXsnROEdTg0/ksBvyxudLymK
oEdhTztHEzhxoiUlW0mE0mvHrOSsFX9M2MFSFuHCheiO8GMEWtot9f6zRoSW
pNB+Vy6B1e//ZDdpiGA/f0h7+QQB5dBL3CItEbqigwIbeQQk5kNszgYRbgcG
jBonEaCqG/rZ6YrgOeBr9jiZgGqbjYO+vggDVDorL3Wl/3BWb2mjCN930hqp
6QRk2fSZJ5tEGLL2Wvj0M4FST4a0Y7MI/wMTC+HS
       "]]},
     Annotation[#, "Charting`Private`Tag$50906#5"]& ], 
    TagBox[
     {RGBColor[0.772079, 0.431554, 0.102387], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVxX081IcDB3B6KTZNaLgSRtaktExYQz7l6MzzSTZy8pDuTK5FQlMOtfOU
h7ieLs93HvJwcR7Lcb49mKkQoluF3zy8rmiu1O2Ufr/fH+/X2ziUST2ySklJ
KfR//n9dHvnfEaV/HH/cmjTiay93XHGp/z6RrAnN0zl1mnrrkctQfyNib8JF
JYGkws4QWi7CUVaEKWSLQyL9vaYIeDJlVcE2x+BJipuTwVZQhWxKmN9OGEJW
1XLXAu8+GJ2KjbCCCZ3SuWhiCccLnsMc9m7IBteZ5MZa4+/R6VVVV+zxfG14
VU2GLdhGSZbtNY74bXn8BZOzG4wxv/vHZvYh4zb9XexFBzQy0kxz1VzxXUT/
0EnuPqhoDNwi5f2IIzF3UhfSnXCwUd+3dIM7Vrjx1q9iyVhWNLIazb0gNx3f
pkN2gcm948xIPhX9WOzl9LniF//i/qzen7Cqf6lO5uyNgq3OpMGsn1HfoVpI
avCGaFkapuMTALqulnKtng80S2yWr48HokHR8v23sz4Qzj3Y2iClwb8xe944
2Ref7M8Q9POh+CHhaeco9yAoD2u/yqyig7U9q2Yl/xAeaxtHLmgz0NvKiW1+
fwg0/8ImnyQGJuqTqw0PBSF2IsmFRI3EEW8bRw1TGkpknlG85V9wZpermC0I
hlx3sbXLk4nusdfuZ2ghGG64+k2i4gRymgXZgX7hOCd/LqovSUCzqbJq6icG
pqZJhaLVqdinNaUuXPoVV8ta98v102HWtijekhGHdt2oXeSWHASNtzG5lNP4
e6QnePJ1PnzaWmL7s87C6LeAnmoXDiYnBPv3qKRg5s+oMEn2Fah7xV5MN0vD
17U3fflcLubesi3TNp2HHl05JHG+GLZpo6ED37Kx96n4+ZuUUoyt5WT162bA
wQ/bg1XL4e/xxLDyVSa095DVHtVWIGTdHtnDqWyYMCcz1W34mHIVao0IcrBW
Q6XIfaISnz+6UJr3KBe3lcO9elKqUZ4Sg21deahyCKpS+ewGMruzEg24+Yie
leS4CWqhPm0lsk6/iOPjbwyD7epRnEAluYQXwG13nv2p2QbcrD6hG+RViPfv
mj9M3hWgeFteitIuDiKWXtcvpN+EnFwQGq96CT01UVVJtEaMiM7dkUgvwV6l
ZUlm2gSLPa+mEzsvg6RWlJ8y2wStHcurTmVcwWRPb1xBhxAPHq+uczh6Fb3G
3fVZZ5uRdCFpZNjqGri2CgtVvxacmCSxtTW5ePx+6MO1La3QHM8tezLDhdmh
eCM9aSsmapQ7XtZdR+Spw+lfitqwQ021LSWtCH02Ovvlqe1wi+4U8N2KkV1C
P/DOuQNxzmTa7JoSOIhXquJIt9BVuznsBrME6udDfJ79dQv3DZbqpgdKcE9Z
de+9m7dhpkud6NpdimLjhLTC+E58cXRUanatFJKX6Xpn7UT4K9LisLWiFCrH
ni1KNbpw/6iDxeDhMrjlfzRyGO6C5XxIi0l3GRxtBX1qft1oYJ6UvN1SjkzF
TgslUTdcg2S/Hs8sx8y6y+1SIzEGDFx23J0rR2tyCUM3Q4zZkRvdGR4VuOSU
42Q9L4aGUrOH3Y0KVAd6OId49qC8+4VW+xoeLB/VRj5s7cGGTduHLjN46C2Q
te/YSOBw30th8l0e3NX/WZnWJ2Cp+OC98T4PA6x5p+sGBIQ/0GaaenkYj557
oG5MYOOd9R9n/uRBSnk+OfcNgd/7aU/dh3j44uMfn5VbE7CY1PbTe8ED9Ujp
z3o+BK54PvOqlPMwKikqekglsHTHQw8KHgKo3P+cO0DA5pbp2NgyD2GOl469
9Sdg90rN5fNPPMRuyGYN0gjEdXQORa3mg/MgvjozisAwPUawU5uPjeS4hX3R
BNZWFmn8sZ6Poo4YKwWTwALdLiJEh4/KyuhOegyBFat+5XwSH22s8EHnRALa
a2hSmSEf9vIQ3Y+nCVBYZl9mfsWHODo4UJhEIHCJYbvZhI++wIBpExYB8+QD
kb5f8+E55G/+NIVAmWwl4eUWPoYofsy8NAJHdTanpprxcbCLKqScJ7CGIz6n
b86HxNr730+/E5g1lrCatvHxX6m99Z0=
       "]]},
     Annotation[#, "Charting`Private`Tag$50906#6"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.885766768691104*^9, 3.885766791478649*^9}, 
   3.885767058849127*^9, 3.886245468658724*^9, 3.886245586103991*^9, 
   3.886245656914962*^9, 3.886245724917495*^9, 3.88624581011654*^9, 
   3.88626617358628*^9, 3.886266205972419*^9, 3.88626683165633*^9, 
   3.886365643258987*^9, 3.886366038184635*^9, 3.886366073592667*^9, 
   3.886366129972883*^9, 3.886366229066216*^9, {3.8863662647884283`*^9, 
   3.8863663477904263`*^9}, 3.8863665933579607`*^9, 3.886368584789682*^9, 
   3.886368692783688*^9, {3.8863688169061203`*^9, 3.886368831047207*^9}, 
   3.886368946686125*^9, 3.886368980880231*^9, {3.886369661777095*^9, 
   3.886369686385244*^9}, 3.886372266324873*^9},
 CellLabel->
  "Out[543]=",ExpressionUUID->"d2aff15a-362d-41e5-b639-ca9190d07350"]
}, Open  ]],

Cell["\<\
The following is the Worster 1990 model for viscosity of DiAn, based on the \
work of Kirkpatrick AmJSci 1974 who actually did pure diopside.\
\>", "Text",
 CellChangeTimes->{{3.886266218606999*^9, 
  3.886266288056203*^9}},ExpressionUUID->"09fcf3ab-1b23-4213-8534-\
9c402e6a4bf5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{"Exp", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"12", "-", 
      RowBox[{"52.5", 
       RowBox[{"(", 
        RowBox[{"1", "-", "x"}], ")"}]}], "+", 
      RowBox[{"62.5", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"1", "-", "x"}], ")"}], "2"]}]}], ")"}], 
    RowBox[{"(", 
     RowBox[{"Log", "[", "10", "]"}], ")"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.886245167868713*^9, 3.886245197543746*^9}, {
   3.88624528766368*^9, 3.886245395391334*^9}, {3.886245578230942*^9, 
   3.886245581267995*^9}, 3.886245650629347*^9, {3.88624580144998*^9, 
   3.886245804683839*^9}, {3.8863656268137836`*^9, 3.8863656296222467`*^9}, {
   3.8863660244302273`*^9, 3.8863660696616497`*^9}, {3.886366104453454*^9, 
   3.886366125832254*^9}, {3.886366224094679*^9, 3.886366341676311*^9}, {
   3.886366473477892*^9, 3.886366504503036*^9}, {3.886366579931032*^9, 
   3.886366588660705*^9}},
 CellLabel->
  "In[544]:=",ExpressionUUID->"b2084bfc-6741-4b3a-a234-020db4a398fa"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"u", "[", "\[Theta]", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "0.5", ",", "1"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0.5", ",", "1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", 
        SuperscriptBox["10", "12"]}], "}"}]}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.886245443282016*^9, 3.886245449680969*^9}, {
   3.886266155653664*^9, 3.886266168654248*^9}, {3.88636853877094*^9, 
   3.886368577054249*^9}, {3.886368640519495*^9, 3.886368686577325*^9}, 
   3.886368805902541*^9, {3.8863689244144077`*^9, 3.886368973836842*^9}},
 CellLabel->
  "In[545]:=",ExpressionUUID->"38fed2e0-16dd-40b9-84b7-1aa1d91d9f07"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVV3k4VV8bRRmSKdGASq4pmYlI3n1OCCEpJQkhU4aSoQGhkjIkQwkNpkII
lTl+KVGpEEKKey/Xna9QidJ3vr/us569zz7rXWu9791ns1eo03EBPj6+omV8
fP//NTx5TpCPbxxWtdrI19cZo2L+X3cw/nE4xTuacLDAGMlknja6KDAOq1Na
769OMUazz4L8hATHwZN5yp52zBjVLrq/ERMdB9h7PJYuYYx0r+xKXS8zDiSh
+kWTECO09a6YrIH6OBS2tl7lom0oV+dqZfiWceAIP9k/o70Nib4QtKrTGIdD
uq5PJDZsQ0zqUqSJ1jhoK2YVpP82RKVbvg+B/jgk3ttp/eyJIVJ9NpBvZzYO
f5cs/nhpGiKl7rvKvo7jkPq0cfG4tgFSsFR6275vHLicsiNzGwzQ2taS0E37
x+GR/L6XueIGSKy6omnIeRyy+v5brcPRR78ymxz3HBmH5b8Vqk5V6qPuI4PR
ur5E/YwZO9DXR5Es8YGF8+NwXaPpjLidHjrpc+PcoZhxEP5O/+NirodOfJVR
fBo7DtGaRvsbdPWQ50f5E8Hx45DgvLBUI6uH9jzZwj+eSNQTd1shbVwXKZ63
1Oq4MQ4Sk2H+PWd10TvRmMvXH45DW4ZOaexpHTT+wElqoHQcxkLs/Te766A5
XD1PrnwcIrsHqePWOmjj+U/VJRXjwH7nFfZsow46ydoy2lwzDnpRQv9Vv9NG
a7oH9enN42DpEJTwUUMbeaXqjqOecTimkWX+7I8mitoidCKxdxx+ePsIljE1
UUrHlx/v+sZB5sURlWdDmujZ38SVhwbGIfs/m1nZZ5pIJOSrUfDIOFg0O/j8
DNFEVQ7XUm9PEPVbfNi6bmor+i05YTozPw4ezfepsnQN5NA846W3QKwPXF1S
G9ZAhb78yScXx+HA2I41Dm81kO3zjV+4f8fhpLfk6ECFBsoPPHyeJUCGbol0
jd4wDYQ6PjRPiJEhPfbMbS8BDZR4rtHssyIZ7p5cE39QZwsaUenyWaNEBorz
wUstSluQTu9gijOJDCLpfu8M1mxBQ2pzo59UyCBMZg46/1VHWwa0Y3o0yICC
v3ce71ZH3TrFz98YkuF+zPI7xcHqaPVkmnmLNRmmkjsKFFvU0ExU8ek1tmTY
02lS9qBaDfWJNpWe3EMGnY/Lj5uXqKF03UlpFQcy0FZ0lxamqSGJ6B1TqfvJ
YK171MDMSw2JrKZf93AnQ9StY26FYmpoCWFkgXAyFK/6fKApWBV9/XRwzdEI
Mpw8WOTq5KOKWnyD9tRHkmHWxbLnt6sqOp9y69mJs2Tgfl98G2mtiuY/c672
x5Jh7oOfhZqyKpoNydV/cJUMYyqbRg58U0H0/NmLNvfIUCTlIxRzVAVRnI7c
TblPhsY2t2N+B1XQF5GXDR8LyKA96Hj/6F4V9CE8g+NcTIaKCreCYEwFPbHT
O+RVRoZyVtpNOVUVFPsnROP8EzJsT6tou/ddGZ2pHbRofUoG95JjjbtYyuiU
v7kHfx0Ztm7w+DMzoYx8+sWzrjSQYZWtiOnpIWVkW1HxN/M5GVw03+T+blNG
sm7MnopOMsC+VcfW3FBGEtJOLF4XGSaicgcDk5WRSFejoMFbMtxReij29rIy
WjC4atrYTYZtHsG5peeU0fhK9eKOXoIPz+TYpI8yqmg+HvXtCxmaH+T9OLhD
GeHy5A1S02QIu0DfTeeR0OuGvZYh38mgua1jzRYWCdkebA3qniGD5DA5JoJG
Qvtv5DUn/SCD9znjn7pfSchH+OBh/kXi+Q6T53rvSOjK7NusGUEKrN1/+sCB
MhISu2HS4ihMAfmPsiZ7SkgoXbuUWiVCgbnz+y3sC0goJ+Cy/omVFPiWNPw0
7DYJlY9BD1WKAo46Lfb4NRLqfvd05YA8BfyNtdbzgknIMYBkYLCBAsZhr5Ky
A0loQCjD9cZGCpwawP6z8iOhb1homf1mCkRvaTv1xpOEePVbdr9WpYCikavD
3gMkFO6cE6KsToHj4juC9uwjofkZoZsJWygg9jZ0i5MDCfFpT0yYa1IgKWwk
5II1CUkX302o16PAjeSyYpedJHQTEy+XNaDAy1dFO8pMSUhu7HzvaUMKUDlt
9su3kxBJ7rCirjEFVBxT+ob1ScgwfXVrqRkFPuidZvNvIaFD0Um/b++mQEjN
uMiutSSksEW39bs1Bc7FTE/hsiQ0PvA5wcaWAuvPtz60W01C/jpqYr/tKGCx
7/7zGEkSOkPp2OjiRAG3d5tnnEVIyOx6MPXxfgpo/T2aWCxE8DeTLRV2psAx
bHHiz3ISSsr20as/RIGcMttfb/kJva2XWaw5SoH6Tcpao4tKyO1HuUiIOwWa
1tF6ziwoIcVCp/cdHsS692nNTb+VUOliwcEILwoUZr8QuPpTCQWV2si/96bA
79T0s3Y/lJCu8/cx5eMU2DZjdX/9nBJqqEIB/X4UGLm95WfPdyUUfYSupRlA
gdymkc3N00oIiaTPXAykwMJ89M8qnhLqOjZ23iCYAiI7sKwnHCWUInEFJYcQ
/rq6u75iKyHHZm1BaigFhjq8Xn5jKaEhmdi0jDAKZLQ1ntFlKqH8Fyr7macp
cEkzpt2PoYQ8Q96vxSMo4H5nRXkpXQkpy0eM3o6kwCsDM+3ZKSVE71Qo+B5F
gcgmPofdBK4If3Xc5iwFzAv2iZTQlNDJzUEaBecocHLR2E2MwPPnmp84xlDA
I0r//a8JJdSi5n2mNJYCA8MF9PMEjusX3ckXR4FPmkdzVxDYIr6W3yWeAiY3
6sj3qUpIRNv19eMECkS0SrRiBO4e4U8WvkSBWP4BbTZFCaVfKdvrcZkCZqZX
DO4T+IDhPpn6RAo8/FLz7iiB15HnhySSiHzReT+VCTyaev+O71UKnP7JfPqD
rITum1p7tV4j8q12RLCHwD5TPNU1KUR+tITpTwisnnWLFZxKgaXdr9wKCMxG
UN2RRoHdlpHeOQSu5tDCN6RTQDpO8s//cXhumknEDaI/x9JViwi8fbfRUncG
BWpGVlDrCPxn9mu7chYFzoTl6PUT+L/7l69EZxP1GwdILxL4kr2WXf9Nwk+b
+stbCb7WC/1SmjkUkN3fnXycwGIPowcu3qZAQwafShmBe/Yr537JpYBVf9O+
nwTO4ut2N8gncORJaXtCr8OVp0nJdyggGFDiV0ng7THB79/dJfIY94QmSei9
xt4vSuw+kW/hqpZzBJ5TOLbZvoACcj5KTA6B+9iu71ILKbByRVlowCRRb8uB
iA9FFOCLmrTnEDgtxWGTZAmxbp54+Szht60mfjr9IQXeh9CWPSLyof5nx4be
UgqM7ze1cyTyJNS9rXNVOQWYyfLf/xK4/cQW+cwKoj9p6gfCiDyaPpJ6mV1D
gV2z+/LaiDyvOy8aPFhLAcvUEupDrhL6abt87dqnFNCM5KTnEPmvZf4KzKmj
wINTxzZlEv2yVWNMOq+Z0NP8gh6H6C/RhaGWLy1EfoK7X6/7ReT3TZ+vQisF
fAuh225eCRUHvG668x8F7EKEtD8Q/atQVul1v4PwXzwj2oiPhBbOPBQjvybq
0drZXkL0/5B1Qd3mLgrE/JRK3LCMhLLoWaJFb4n8h69MVyXmh7h6dG3JRwr0
aTfR7oqREPtXhButhwICq98cdJIgobedoUJqfRR4x/niLCFFQol+3q6l/RRI
iHQWvU/Mq6UHtgKPhon+sV6briFPQqORFhWsEQrgb8OblTaQUJOV+UHNUQoE
ZmdeJW0ioUiaXnnlNwq8Wde7DkjEvFZZ71RNpcAaVYF937YS87B4quAZh5iX
J563tJmTkI0b3WUVjwKTOzeRzmEkVLuaIRk8TQHh0E1McwsSupTAjFGepcDz
rrhHUzYkpO7JOZw5T/BRUDxX50xCJ+VnV51aRoV9r71O2YeS0HDfbFf3cips
rz4vXR1GQvi1uQvqQlQwbBHw2hBJQjK/f3DGRKggx9NUWhdDQg2D828cJKhA
yUh+mZZMQvwZS/Ga66lgOfSt7HI5CQXa/NueJEeF7OVOk9lVJNTPxzdNladC
pe4t7ZpaEnoQyu+et5EK8iPc5Suaif9X++WmospU2LTfPsa8m4QyRERnprSp
sPFNaevraRJSuiDrVWRBhV82jSnxuDIqOey/NGVJhfKj7Nai3cpI3bA5V3M3
Fe6G1gX22CkjbYbnp2c2VMjYmWRnc0gZ7dhfYdHlQIWtkQuiTUHK6IDaLjWO
CxWmjZ8b372tjBI/nmQbBVNhs8QmdvmiMlpR/jLpfAgVImXj+7YKqKCUS2tU
/gulgkJ5oHmDiArKMGk5ahNGBZOgH92zsiroXpHwxyNRVEivYXz6qauCGqLu
1sTFU0Gl9B1sC1BBjE3dkW+zifpEi22/UlRQxQqXapFbVFjY5oo3MVVQyCyV
YZVDBbHfGdlFMypo7vWi28tcKnw6v3A+T0AV8YVs3dVyjwqvdOSSN5NU0ZqW
a5JVZVR4XOUjesNXFeEuNqU3nlPhoLR79YFfqkgQHyB/bKXCiobDpf18aqhr
6zF5if8I3KiQ7iWqhuz/nUm71k6Fdb2n/Eo3qCGXB6XhFzupMHvwR8dXCzUU
OieMRfRSIUvEYd2fbDWUn/562GWSCkIGn/lO71JHzGOHFWtpVPi6fELg9V51
tN2A7buSTvj5SKRlq5s6+jSwau45kwqfN9RU6ESoI1H5oxKkaSocq/73T75U
HZ0pmcU5C8Tzz1BT7qotaH+TYkW81ATg/9Qff5ndgu4nP5kZWTUBWg6dw0+I
+zvXzcrEcPUE+M+jwYJVGujavxOvabITUEg3eNqmo4H+s6gf3yM/AZ//CckL
BWsg7Y/2smtUJyBgIGF2kKWBRCfOxZbvmIDyr2sb3Be3orAhj3K22QT0rN1o
8XulJhrpthjUNp+A+78WtcoVNFHZMwntJ2gCdmdop+w110TWVwq/NltOgN+d
vZmPEjTRFY23Zu/3TsBQfsqpdkkttPzU+kWe9wRwqubGcTNtFHR8SVX/+AQc
/WkYS92rjfoPU53CfSfgt8lxz7ve2qgIryyf9yf4WFyP9U7WRrgM5sofMgEa
D+8sao9qo7h6/0bpMxOgyO1KYibooL9/66O2pUzAoz1Lt5bCdNH0cFF1RuoE
WH0MlL5xQRdRnl1n8NImQOFxQs+2FF3UFex3pPzGBCTEjV2qK9FFGV/XwKZb
E9DfFzZwa1gXqT+PEFpRSOjJ1/ZrxkIPHYg2yB6tn4DkKifXrxr6aLfLpg8m
jRPw1sNb4ruJPjI1XCl8q2kCso9JDsrZ6KNNbMqZfc8n4FVWzJYaP33EPJLh
9rp9AhIrUh6zS/RR3I7vpOr3E1A9bp/YqmqAHi1U1V6kTsDTrZvpA6aG6Nvk
p6SRiQmQ8ZE+c9fOEEn1zrvr0SbgSGMlN87dEIU/xFeO0SeAPv/MICvBEJkf
GPQ25RJ+lp8rnX9niPoq/8p8n5+Am8MuOd+8t6F5T7sod4lJqMdGTV3LjJCG
3Sn7Z5KTsK+hivq01QgdNb5JEls1CRbrhrR0+o1Qu/j4x8bVkxA++zH84j8j
lNJ0Wl12/SQ4aKy7Wn3IGG2SyRt+R5qE/WpKM3sktiOrTsYOE5NJOOy/N3fN
TROUMeD564zpJBj0kgRlqkzQGPVzbcOOSRhV1Nph8NoEneF7vWW7+SQwDBdG
P/80QeUmhWuNd03C84frhtsOmyKJCtdZQ/tJqBVP49uhugMNpr8r1z02CeXK
PWERX8yQ0j3c96TXJLS47jguN2+GQisbN1d7T4KTQHUiXWYnEnn7MEfHdxKO
Bxl78u/diUyXXUrUPjEJ/858e6vZsRPdjTDz0oyYhG2LSnwPG83RcdfK9epX
JyEz7i450wshM8etlIRrxPt+tWMGEQjJWJWVf00m9DET2j99BaF2vZIdmWmT
8JK10v1VJUKKK+4cXcqcBLPC+/wjCwiN1qcU9N+dhC3RG3azczHkJBu0Jf7p
JGxftkxkzxyONFYyv395NglG7lsNpVbuQvz8/k1G9ZPwZrv97YXNu1A123sP
u3ESPtctKiLHXUji1ZHgQ22T8OfRhdfZVbvQm7A9NVpvJ6EqN8euPMwCoR4N
k5GxSVjuLM1Nk7NC+owF2gbyJAzccImpMLRCJIF3WccokzB4vJEz42CFhAxP
TNMnJsFd/ELG4kUr1H3r0cNfjEm480xhetm0FXLx0FwjOzcJH0qCF/R7dqOT
XK3ZvSI0mM0Ze2nyxAZ5CS/dz1xBg78OuR6in2zQfsUPDp9FafDkU5uQxKwN
MnIKKXcXp8GD9zH0xwa26M+zqmPB0jTi+/axqX+9LUqK0em5pkCDJd+hG/+6
9qB7YnqVHbo0WO/8ZTp7rQNSEBfxV9anQcPVzOj3mAPKER9TumhAg6SLIXsd
gxxQukRqDjKiQeHDT6qm7Q4oQYpxsXkHDcyTK048PrUX+coUuNZY0SBTVSpv
INQRaSqsErnjRgOaS+g6y6V9qFSB3r54lAbGK0X7nq5zQiob2mJcPWggtRz+
HjBwQhs3Bs+u9aLB0LSKsXqAE5JSfPMtw48Ghws27lP57IRmSfHPksJoMH5K
dnSgZT9q2vrdKzyJBuGfqD8Va51R5A7+5u6rNDgxqckY6HNG+ntWrVZJpkFE
afP+zllnVB6o93IglQYBeQ6Jx4wOoryyUyTjTBo4bRB1jG89iKLVZqi/79Dg
2/6IKvLIIbTdWMDM6R4NEiJOb837dwjNWUlnld+nwedjkmdvK7ugIF99iyNF
NBjtqnjtG+qCjpaEFbeU0uB90EvjAyKHkTlp1ufCExpwh78b/Nvtihb0lz3/
/JQG8yKbVtidckV1+GpZ3ToaIJFeu6lcV6TtZdAx3kCDm72PG1byXNGmgtMq
eCsNBE2NSzbmHUF8m+Yml7+hAefZLgGhZUfRC7kfvte+0eBpyeLFmJWeKMf+
UPG7MRr8Fshs+knyRCFxjWQxMg2y23ZF5Jt5IjlajNt1Kg0kRdXHQkI8UViN
kFMWnQYt9otXTfo9kbL1up13Z2hwT+fVssbSY+hKxI7VtUJTkNA/jJPOeCP3
0juOs8JToEe+Jc666Y0Mv/xLNVwxBZLrju5gP/NG4+iVSP3KKbjWFJabP+eN
tovZLzVLTYHiDv/c9nAfxCx0Z3TITYEpm/UqPe44su+JbxvWngJtSXLFyBs/
VMVkrjDSnYKLghq2MnN+SErwwIEMvSmwjNnmV7LRH/WZqDJsDaeg1wvHe8P9
0cGit9ItJlMwcFlu4LhKAHKPXO1712IKhLf39ezICURt6dHVC5ZT4JgmGre3
KxApPppcOLh7Cqat+TZ3zQciylj9dUlbgo9O4RjL9QTys3FrjNs7BXjTkM9a
5SB0UqFEzNt1Cr5Gfhb4/DYY9RpJHGo7Qpx/+sXeCv4QpL8vqkD+6BREffqA
/dweguYu2xj1e0yBPqfMlV0ags7yOB6Wx6fA59qd0QtpoSih3eiJ2skpeN9Y
de9+w0mU5d91mH15Crza7YeGP4QhyQcXr+pdmYKH02++7p0OQ9eo5o2RSVNw
vlB+laj0aRTv/mwdfzKxXqfsbX/wNAo5UPhZNn0Kbt6prVWknEa7UbQz5E6B
xGLax0ThCNQeY3z5Ut4UhC6zMLmmFYHMm2eevsmfgtTo34Kj+yOQoVGAzIF7
UzCymZsSXBCBFDUPfvIvnoLVGy0k3VAkWliruy+jivBv87F00bQoFOHMivv8
eAralu3utG2MQtMZD6oVaqbgVDvLfHIiCk2Jb5B6+GQKtpc8VdfceQb1LxP9
2NwwBbs2yo2tmzmDKnlUu8n2KfhOXXTWDDmHPDtvWW//PAWan77yCp7FIHmZ
FNmoIUK/2r15HewYNOgZR3k2TGB37D835Vhkt+AfbTBK+L+zPZuRFYuMtUyr
tclTcO4cfTD3/AUkkTm6VoU1BYcOXG1Xq45Db8Z6JrzZU5B/XW8q+2McuqTZ
UVPImQKdubv3rHhxaKGjco/i9BTYntfXD9OJR5PzsRfkf0yB3bXIJYEn8ajF
ffPUqn9T0Fc3i253JqCoR7JPHfnoML5N4d15VgLSn18Rf52fDvzZu40eSV5E
pTdm5MSX0yEhXTt61uUiynz10kFkBR16Nwr+t5N3EQVq+NYtSdNhfo8zzX3r
ZbTuZ/llliodlrobf/hOXUHuGg88LNXpkORZcum4XBIqdi8wubeFDm7cdwe+
2ichnc5bHCdNOpy6yczbWZeErHIuOTfp0cFGtKiBmnIVRZi6q1w1o8MHKycr
hlMyago5/I+ykw6PGtQaZdOSEV/RgWEzoIOZUUxmz5tklLJyT+o0RoeH1Rl7
fuEpqHjU+Meh3XQ4ubfUSt8sFfXFrupQdaLDj9eiCg17rqN1T8Tuxe2nw3Bh
Ynl8+nXkPiV8duQAHZKNXdtmBq4jxt4lzbRDBJYp0g3QS0d8m1lZP9zo4MVI
C3ZrT0faL1/5vPIj+L6v4bj/uIGShc8s94qhA24wrvPxQRZiune8LI6lw8uZ
yhrfvixkWyd9ceoCHUj8O3O8l7KQ6PFK/uAEOgxu4qTcOJSNrrVTlqKu0OFq
5ItzIxI30dUY+9+pN+jQOP+1WCzjFqL35zb0ZNBhbzjDKerVLWStSY9anUUH
SvDyVJdft5DIl4SfOTcJvyJ3HbzonoOStjfMFuXRoST10Ks8/dvoyqwSt7GE
Dnopq6V4M7mIZnuy8s8DOoDahZ/OWnnIqvB5EJTSQeiuWMVu/zwktM+F9aqc
DsZ710v1juWhxKpU+sfHdOh8c8KzaCAfXQ6Yp0w2EvrlNmxM/3oXXRx7/3nV
BzpsZIbpfHlagEI2KhWHfqSDGDd5+NpYAXI5GnnyfQ8dMq/HUctFC5H2l40r
rn6iwwsBFXHMqxCNDIaa8g3TIbokdhdjTRHS/bDqDpdK8JtOd3W6VYzkxX0D
7CaJ93HW9xW/KUZCdk3byml0aHYwsI75U4y+vPH6cJxB+N0prPHauwQldjz5
N8qlQ8tLzi3K9gfoa4uz17vfRF6ddLdXLD1EXYtl2lsW6WA1eCpD3rQU1Zou
LST+oUMGz3pRI7IUJTU8yMT/0cEBr9PWmy5Fhk9/vWpczgB3q4V3EuwydO3R
bbVSSQY0ki+bti2vQMa5X9mXVBiw8Pjaeg+fauRFDzw4qsqAJ9YOOXtSq1GK
0XybgToDqlJ96x7UVSNyn3QmRYMBNxXlPlSK1qDkldYmSJcBKRtxXVJDDRqL
rr28aMoA72+1MwNbniDRd8BzMmOA0DoPS1vXJ8hw/XuX8p0McHTqXsKTn6Cr
z6a2HkYMMLT9WNnJfYIMuAp9dZYMGL36YrtNy1N0xePKxjBHBsReQz03w+pQ
baVM0pt9DLi8pjpUqqIOjS4WfFfczwA3kzXv1Wl1SC+n5dVHZwYYN9RI7HSr
R196vgdqHWGA0wtXIXXHBqSDu9XTjzMg0U51eatLE/qsrOfocZ4BOk8uLLIX
W5HCTal796IZIPUvmduyuw15Ck9zxmIYEPg5w+pAVhtiMqqSPeMYcNB9jDO9
9j/0p0qz89hlBvgOm1840vAfUjRRN/e5zoCnNM0NZeLtyKdcOLU4nQHmffxv
nE3bUZn81JeJGwzo8FyUO+vXjvSXSs4ez2KA5ilFtvrLdmT5kvTM9zYDbtF7
G5ziXqJA+02aAUUE35nKWeHVHaiqdelcWTEDriW/MHy1uwPN6nx7wyhhwEz4
h3HxmA4ULX3HL7CUAcssYl7uZXSgtM9yRScqGTBicaStqes1enJsjVxIPQM4
/7ZLPSjuQvN9P/yrGgj+QSl1/yhdyMxioJ7bSPBL5WkKKL1Br1WznENbGGBU
L6ewovANGmKtunHyBQPmtXMtnMvfoj8R4iKn3zEgatDwQMXXbtT9Ydb7QzcD
8hIuUTu13qM8tZG2LR8Y8EJFVag39j0yHX4QNd7DgA/b/R7Xkz6gM2Zoyn6Q
AZHt9501zn5EcwKnX6uTGYBgdE2tRi96deSw0iUKA2p5yl2Trr0o8ynEjlEZ
4Ow+O1mW3Iv0fcWMbtIY0Ar215ncXnTyTUnJMjYD2oft+eNb+hD7+tClbz8Z
kKZ3QbU2oh+10FvHTeYZsJmdqNzzuB8lYyVm2b8ZEP/l8Mq7zH6kMXtqbs8f
BuQfKfNw8hpA/gdX+jTyM+Fw/ikXN9dBNKFgvitLjAlFNMcSkdAh9DRC+R5P
nAm8gqBlO58NoYsfRBdtJZmwtJsuZP1nCCnFfa7ll2bCS5uUT+yUYeRJPakU
upYJ+lvPurc1jKDRsiJ+WyUmYF6Gj9bs/Ipeq4ziaSQmCCmXgkfKV/S4QOZS
nzIT+C0ydqeMfkUJty8LHVFjQtrb+yl3L3xD6lf9VgZpMmFVdvP77z1jKMJ/
q2yaERPmpHx+hu4jI/cJ74N9xkwYjWL/oMaRkbVn/q01JkxoM8J/kqvJSMFF
XO7eDiaI1hjcTJCmoPbdvI3ViKhv/HXn5VEKklB7ot5ny4R3yr59m1Mn0Hwh
K2CNHROyBn6VCL+eQOSNyo9c7ZlgzSGdyfo3gZ6uydai7mXC8pIVGo/CJ9ER
oSj9uQNM2O7s9VjIj4ZKJ03N1ngwYSwieYVbCB1lHDsd4+rJhJyy9LGaWjo6
//VR691jTPgOuHzfLzpyGFDA1HyYsGnY+v75Sww092rJcnsAE9xk/SsbHzIR
Vty+1/U0E/7myH0xWMdBgTqG2jnhTDjuWXuoy5+DsppKxAYjmOBi3kixbuIg
Wg8xPs4wobx2Pn69Jxel/HXYZRPDhMg7287a1vPQ4MEvRtuvEHzItZV79nxH
/8h2spFJTPDXn6gtP/8dqQe3zj65yoS7Jl2Pciq+o/MXC6p1UpjQPigXoC41
gxSr/TTUbjBh4G0J2frbDApcMbdhTR4TbKi37uvkzKHMrON/9uczIVFxtYpr
3xxq2fR55MYdJlxDkqv0xX8gSaOmW+L3mXAlpC1v5tIP9NQrbpVgCRM+8sX7
zUb/REvNYoJzVUzINC5aOZs+j9R2x1L1qol6/k6d5/88jxz7eC9Cawg/PKdu
Ujb8RkX0vljWEyZ4bdXUU6v6jaxlb89TGgi+oUVZy4cXUEaICru3nQmpfdRT
Ca5/kRX5zmeNV0x48Ni5AK/7ixb2r3l5sYMJD7GUGd3VS8jbRDh3Wxeh/9Ip
29KeJWS4nLH79ntCv07Duzu0+LCB2xVFnkPE+94eUm8a5Meuialcbxwm8mmS
OSwpKYCZX7hzTvoLE74cDnmtaC2APfBJ2/fqK5FXQYPLsS0CWKTOyX9qVCaU
6Xn9C6tahq3rMDjC4xB5PXRxy6NGQazbuMLKmseEDsGnoQtLglhcubJ+wTQT
lpVP6epbCmGM67Ir9s8Sei6x6H79QliT66+6unkm4NL4ZuslYezIdJN07DIW
FO57VEZNEMWkvA2WBpezoML2HtlpSBR7NfCIoSPEgjeMARuq9kpMszn/P7II
C1Zw9z6OGluJ/bkcG2IpwQJZXVlNioM4dlcefyu2ngXeFx8sWG6TwoJqxPIc
5FjgWL9d7KyvFGa6+/OJdHkW8KLn+61zpLDPYUHiMhtZUBb17+nrv1KY1Nub
jvIkFliGy7lF9qzCEs6wBrdosUCTukEq8+ZqzFGi7uEJbRY8s6/Kz+5bjW0s
jjtTqcOC/vmHTBdJGaz5o6ycnj4LuJvM6VuTZbA5VXR0uzELKp1r1O1SZDHf
gSyqFcYCBznBLSN1azHDEx5Pk3AWPKTQ8a186zABfo3Lb3exwC18uC/Adh12
T7NN1cGKBVF6l6jN5HXY0EVGgPMeFpgyvshLKMhhtvrm097OLNCan7nn3a6A
resSeVFykAXK3y7nKShswGhHP92YOsSCe6eTby1EbsAuXvU3OOHKgvSBGa1t
uhuxlvGMyDAPFmy5fW2LbM0mTOf61J+4ABYMy2p4rj+zGUsSaxz5GsgCG/bE
pdYnmzHK1WsNpkEsWB7rUrF9ejOWdUkrfC6EBVbGp4IaTihh82dPs/zCWWDc
LLsNgklYu8+/YYcLLOhIiK/yuKOCyVN66h/FsSAvLtuWOqWChXsUZosksGBM
SMn9tIEqpuZqua/9EguaZN5Ud71XxVL2Jndtu0bU65X/tlFSHXM2XVuvkM2C
Ax5njp3t1sCq6ulZZ2+yoJlcNty6ZSsmsq0pbPAWC1KjFPcdStqKNekc1U7P
JfJk7Lov20YT26hSVLLsPgvi7HojB0e1MLqkThaznAV/bA8H5JfoYlgqX9ju
ChaoBW29v35cF8sV7dtbXEnUtyvVMlheD7MXjFjpXk3kDfcvK87Sw2oXmhJ6
n7HgTnjCzuob+tj5SatTDf8RfjhqjNk0GmLcy4GGSy9YMC19Tpr61xDzUk37
teslC+yztAUDdm3DrP0GYj52EP6+0w637t2GrWF4X6O9ZcHUtCX2bN4Iq2bH
FckMsODbabPa0LMmGCm12Nd1kAXqOfNBsR9MsJtaXVvuf2aBITlW4pmyKRYT
Ilm9dYQFgyWpL8v7TTHb6Tst+Bjh1wwtR3eXGTYx2zRwksGCiFJvjwJ3wA5l
fcupY7JA78u3h1ZdgL01FHD7w2LBtU9Geo/XI6wmwoZyhcuCjQ9iVf4mIiz2
12fu3VkWOMUHi7UGYti6xTnh938JvzI+rj50eBeWnLfunfQ/FpQEFCVJ5u3C
lnaYpbnwscHTnYVVf92FTUZflJkQYEN05ezn+uMWWO3Sqs2LwmxwO3O3DuIt
MTsBHdMt0mxY/fWKowdjN/avJjnXfzUbeMOXGspMrLHaY/SFhzJsOK/7I1Is
2Rpb96KgWXUtG3LiKF809GywiQsyO5UV2DB+8tyHncm2WPTf37BJlQ383pdM
L52zx3QqD953V2PDJaNq8X3f7DGK2xO+u+psUI0VXore5YBZtwS9UNhK7D9d
iAtJ78Vkzo/hcrpsGPwREGLs6Yg9mn9lKWvKhuku94UqAyfMvXTzgwM72KD1
uCPm4yknTMolVijLjA35l2Z3dFc7YZH1xp3SwAb/6Q5pP739GB5Zbi1lwYYm
Cd+omJ0HsJHZ63tWOrBBzjvcXO/SQSy1iP3Idi8b2kbLfR/1HsTQfhuxa45s
0PnxyMJU8RD2oFbgvch+NiS/iFd/8N8hLOxUhIOQCxsmHtxS0JA8jInyXPfx
e7GhkFJm85R2BLNObG094s2G5avmN7TudcMSNyhp1vmwoeF5hYZXkxsmYMcQ
OuHHBnkz86mMrKPY79Kolv4gNggEXm51OuqBGaEvW3RC2GA51XqK8ckDC/9s
futqKBs4lw0z1VU9sWlBoTDzMDasX3C/4VLgidG9stRKo9igdrTthkXRMUx1
4VeWwFmiHpaHnzntGOZz44jA0XNsmDtL/OlreGHjbUpfV8WwAb/IfRf4zAv7
rFCTcT6BDSMj8HB4yBt7Pfj+r0MqGyY7X06quvhiy0P0TpSlsaH8T1h3U6Uv
hgtmDy1LZ0PpkSMOhcv9sFYDtycNGWxo3dY5H1/nhz1LZ/or5bBhw4Oq/s1b
A7AZNYfB6Nts+FIzSeMlB2C6bTW7PueywXtzcSlwA7AKzpmNKXfYEOLduSGm
PhArthXu/1HIhisJtiOko0EYhRyIORazIeaA4cvat0GY4tkPVeUlbDj93Gny
l0kwlv8w+6pHKRvsbqkZv9gQgmUtV4Y3lWzo89EP3rAYil1sReX5DYSeN3cF
0VeEYSNnfv4Za2QT94GemdN7wjB9g4q9pGY2mH2pm/uRGoZRHq79UfacyMcc
M4AtexrDb3Chvp0NCV/OHpfUD8dy9xRnLLxkw6hyfGXT+XBsRtB1cmcHG6Zs
m/miX4djhec6rr3qJPxd9qx9n0cEJuCTP9DbTej564T9q7uRmOtGJ3XZD2yg
FWvI5n6PxGqHhM+7fGRD1rGMojLLKMzL/rTiWC8bROyr7K9+j8LajWxPsAeJ
fi4rrBpwOYslrJjnEyYTedxxbcIsMBobeVm535ZC+OccTF//KRrTj/V+kEpl
g5QYS/DKzhiMMvNhjwyN6JcaUTONtbEYPlpyczOLDcP4kTGrbxew3JtHGD5s
Iq+bzLZs1YjDZhxXmZVy2MT36Gm2jG8cVtgRTdaeZsPzS7dT88fjMIHH+zXN
fhD5W8T3ksfjMdeAFRcu/GRDxp0/ou6kBKyW1Nbb/ouYF03e1iO+CZhXjkaU
zQLBZ/GJKnk6AWuP528/+I8NvXfp6/NWXcLiDzw+dGoFBxynVO+nXk3EfH6e
qP8oygHWZMTRzZOJmFWO+lptMQ40bJyN2IZfwcRGCwaZEhzozH39dp4vCcvx
znL2keHAyCoWt/TGVey84L5n7bIcyBeMqlr8fRVzfyguu3ktB1KWEucFfa5h
yqzE/q/rOTC6oX6dx85krDrs7P5Dmzhwc6/d8T+CqVimjNGTZ4ocUJvViKBG
p2KRdTPSMkoc0MrOCHjxIxUzWzjR16PMAdt3OW6XOGlYZ7z7PlsNDqiytFok
z6dj39J37d25jQNfGpYp/m3MwF7o8z3OM+KAiOljN2fhTKy4v0ViwZgDCsbT
pr+cM7HAdUYf6kw5kFUnmaL9KxP7cU/dXhdxoKIv87ekbTY2hE1WpGEc2Jjq
K3W4MBtrohaIcXAObPo4/ln+TzYWpybfXWZJ7D9y9rr605uY2GPxPaQ9HChX
meySN87BeI5vy+PtOLA78Mdrbl4O1jeTKDpuz4HX5jtawgRuYzlGfG/zHTlw
sFIu/VH/bYzUOmO95iAH/hP2TmhKy8NM33+2WnGMAwMZLRUZnvcwt9j3u455
ceA863bdEPkedkHnJWr0JvS5K0Ar0L+PvbpRZRrgywHpc3/dDjXexxwOXtZ+
c4IDmR/4dML7C7CTwue3bg7mwAvzNH2DjYVYZsNJ9bMhHBD288/UDCjEhuTc
lLac4kBtZerKCYEizHtMf83VSA5gIg/TBKyKscR09dXkKA5YSBRk6OYWY6XY
RimTsxyYM/vW+ppXjHGLVogyznNgPE/khMW9EuyM//hf63jCX5VgSxMp4hb/
PXVSJIUD6TE3vcVROVZVeJHimcqByxKuspOV5Vjv/rNjDWkcSNO22HpzwyNs
TZ3PsP8NDrQq7ggoEajACs7teN91kwMd/H7PfScqsVdb9d4q5hB8tB6q/DpW
hdFGVTvP3OYALbuo/AG5CtME6Rfq+YR/MvbMpKnHWN0yxtOkAg4oNzRHGgnX
YENPv9WMF3JAHntOXba3Bls83l+1vZgD32v0lq3IqcFQV1sp/QEH3DP+qPdo
1WJvU27mW1cQ+v5s2DXt/wQbk7W4JFLPAb9EP9MXSnWYAjPZ16aBA6lMq68n
z9VhLq2frK81Evkb3M69/KkO6/P1FhdrIfpDIbPv19V67FV9wk3JFxwo+GOU
UCzSiPGnvD3j2M6BIfGbb9SDGjFzT+kjN15yYJpzbbVFbyNWL1K4afVrIg/B
tbGu95uw8sMvSte844CBGvtFl1MLRtNakXyomwMwujHo9csWjCSwLzjnPQeK
P11rzTR6juWXj+vJ9XBgNqn0X/PmVuz64r+mDQMccFm86JEu/B8Wfsf8vfIY
8X7+zoEPqu1Y9anEx8fHOSC5q62yKbQd41h+uPGAzIEw2ZnnnxrbMV/u0YPq
ExyYt19xvWDfS+wwxIxtZXDgWovVX7P0V9jN1R0vgphEvtbbSZCor7C+KTHi
ekmcF/L6sadxB2aXnu+vwyX6J+Xk9d/UDsyc3Pxdf5YDz9/tadjt0Imdf7as
//QcB1a8k8y/XtGJNVzdU/f0Bwc2F5697CHWhenpfzlnNE/Mp1b9kpqeLkz5
4sIy078c0BUdkvka9BZboWKyBhPmwq7sFzYlv99jPeNXBO+IcCGkfn3e5tgP
WE7+4Nz8Ci4sU7TO+bX8I6YuE/6pWowLTVZV1lyxHsxq2eN0RWkuZORt/jeh
24tJtP29EL2aC4c7mqbUY3qxgXN2oUMyXLgRs9an620v5jPDsE9fy4UN14Ju
GgT2YQlk5ZV8G7iQ7mGyPrzpE2Z75/TCkY1c2P6yYaW9dD+26nA7o34TFxLP
pLyqP9GP3e/x6ApV4sIt+ftynqQBrK0t9/KYGhc2ijjofXowiP25I8XXps+F
oNuPD0UJjmCvDnvw5Ay50L5352e9cyNYsmzVt8htXJh7NDEn+H0Ek0vd81x7
OxdUWkPGlCa/YNujE8/d3cmFqKPveI70rxjf9oGABXMu+NJ+5s+HfMM6Z0mH
DyIu9NDUiyTmv2EHT7wwltjFBZrjqn/DIuNYhOufuRhrLmRbbpCX0ydjZmv2
TAzbEPxnfFs/XCBjy/tuf9q2hwvCi+EdrA9kLMvGuJZjz4XlZ+XHz56iYLXb
w0KP7ueCVEPj47IuKsZdQ2eYuXNhcmbVI1oHDev5kdQa78GFsI0KA3dVp7Ca
/i2Zrz25cLvqetCfpCksPOOEmaM3F0yPPfbtcKJji+K8NG9/LvyRL8n0WWRg
o+zr3qUBXOj8Ird/QyATa32nu50TSPANV5BqGmFiCVfDyJHBXPjVsz7hxQsW
JiL40+BaGBecffeRlx5yMCb1psjH08T5AS2k96pcrLvd+OvqCC4c6go/2/mQ
i12PO5t4J4oL/35qp/pV87A1f/4MVUdzwf0/pzcN3GmMNLss9nMiF0q2lXSt
cZkl9Ch2UkjigoCNQO6RR7MYrdpS7dhVLujZCzEk+eaw8pArPcxkLpQbnBIn
V85hukxR0lI6Ub8d7wxF/icm/ebRLzyDC1/OK0V5nP+JzT20676SSeRZKdX4
4uhPrP54WsSqm1zglbIkvhb9IvpnVZdyHhfO7y43f7r7N6b4X21+QD4XHg85
pKBnvzH+e/tPVd3hgl2nzPEnygvYK7dsue33ubC2TCAodMUitmdoXfCeEi4U
l7SUm7H+YC4fN0mHPebC/s7MC3YpfPi5578vVldzobWIHSr+iQ/Pe/TpB7eG
C6fPBDsvl+fHv125MnziKRcuqYV1l1Xx496Id/94Ixc+pqReOzUpgF/SfiNd
3MSFdR3s/AtGy/AHCkWXKM3E/lCVr9+TluH0+YP+Hq1caKsfElLQX46H1LTq
HH7JBfU8NXznTUE8/V5Owa1XXIi5wrJQmhXEa1LDVg92EOtGu57f3yeEzwWo
/nLq4gI/dkDoymph/KxSWqvdey4oJFgqv3skgl/KPGqPPnPBIjkibvs2MfxB
vHFb7BDh54Z3Ct6PxfDO0FV6z4eJfhXXjn+0VRwXteuQMR0l8i8totuiIYFf
X641akDmwgX+xMiFzVJ49YyQQxiFC26ek/rS/lJ43/h4WzWVyJ/t2qZjj6Vw
medZRVo0LswfqgxoxVbhuRF/AtVYXDiSYPheMlIab/YeHD3OJnDR49rE19L4
6L5qh2IOF3QYsvk+61bjm7R99DdPE/mvc9EK/W81Xkx791vuBxeKsN8fzirJ
4o8P5SVJ/CP8rR35VGi4DnfTihb25+MBb4HrzHi4DhdZdjTxBT8Puprqq45t
WI97Vm28FL6cBy6P3/evlpTDpQULL4ys4MERswfs5xIKeOtIwpLBSh5Yap8r
681SwAOrvWNSxXigwaUJRGzYgL88onIeSfIAna1euGO0EQ+vLY18IMODGzrk
d78sFHHFpKs//snywF+UkbaQoYh3Hw0MP7yWB2Gn/ASjKYq4ygrNMDE5HojG
/NNYn7gZH/R4HBy2iQea+Xcab40r4dvF63zMNXhA82N+W8FUwScoN6m3tvJA
8jA9qX+/Kp7eEOX1XZMHIUtjaWVtqjjD28SzWIcHoRo3xPLvqOG5TS1HRLfx
wCvJWa/n9BbcKv3OiLcRDwxvdzSt523BZ47HHn5uzINCttdCR5AGvmcVOnTS
lAdxhg2bokK34n/8XjoNAg+e9X5NuJOphZfuLO7VwXgARbn5biRt/MDqy45X
cR4cNv0rSXumjVe2WjmYWRL6MAWn06k6uLvsW5tCW+K8rD9PwvT08BftPSjE
mdATYpfa/xngSc3LSqoP8mAwnXtM4Igh7vDUaMXsIR7Y/ifjzKo3xEdK8nuj
XHnAEqU0DEduw+eS/L3iPXjQUYQqJEWN8ab4vNcvPXkw+UXoRfxpYzzu3AcN
IS8eRCfmriZ/M8bFgwxnr/kQ+vtmB6Y0b8dV9/JdygrgwZPzyg8LU01x9m4D
+mAgsd9k0OKa8A78CfK1Wx/Eg7tj4h8rL+7AkX63zN0QHuzoqnnwKcEMPyyb
U/LwNA+ivHV/9d0xx1O+6HQ2xfDgi/TJfddLMdyp32vr31geKOyqN9i8gOHr
32dfhzge3KqKETW0x/EHrYuHXiXwgD+86PHvXzjeVtBJf3+F0HO/QfSQpwWe
mLtgJ3WVB+bNNfUWbRa4XaZWjdM1Hrz/cHssdJMlPnQp8+znFB4M2BQNKk9a
4jN+HqLkG0T+NC73/EvcjTd4ZoSQMnnwIry802VuN37hcEff8SyivtXe/Q3e
1rjYnq15zJs8UMTvlzyztsFVtH9tncvjQcvDRuMXmntwl7k0e+EHRH2cKssO
37346w+4buhDQi/vl9PBgo64YdlP6c+lPIjI2iOmc9gRl3R3H3rwiNDLrkhu
XnAf/rpTy9uqhgeO5wSU6yOccMNCsmVVLdFvQ29k4wed8MLobPU1T3mQZOej
t89kPx6r95czWceDmwflDukKHsC35XWfudzCA+1GS32deme8MCLuCOc5Dzb+
kK4K3nIQl3I0NHdu44GaCu3FqrsHcfby/GUq7Twok4DotrRDeHHIibRXnUT/
SSaqGlw/jEvjosXLPhH6VSgEZqx3x+MUWhNP9PPA52Qr92+DO879eSrg0wCR
v2h3484jHvibR8PaRUM8OHtFauibsyceJ1vWiH/jwVzti4nK5mM4l+eWXz7G
g3BuniJzvRfu9lbqgjSZB7HhWr4q57zw7RfOWFCoPGjyWgjtQN44j777YzyD
B/eji8Wnp3xwt5eLNQwmD/4FGif5HTqOv73zOGsfmwcSX4NCdr45jj9wWuu6
mccDwcqKfSZPfXH3FtrEf3M8+FmRPTBb6o9338ztVP/Jg16dYsUxzQDc9JRD
efovHox/4fjmPgnAZVXrQj0XiDxd8fjY9iYQT+ALdOpa5IGdgYiV7eET+PTI
hm26f3kQY/rDZ5J1And/1rs2Z4mYH0m/zAvjg/Du65cX/v3jQfvmmt0ZcsH4
/wBeoHH/
       "]]},
     Annotation[#, "Charting`Private`Tag$50957#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0.5, 1.},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0.5, 1}, {1, 1000000000000}},
  PlotRangeClipping->True,
  PlotRangePadding->{{0, 0}, {0, 0}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.886245468714591*^9, 3.886245586141309*^9, 3.886245656949854*^9, 
   3.886245724948724*^9, 3.886245816201882*^9, 3.8862661736475573`*^9, 
   3.886266206022924*^9, 3.886266831705504*^9, 3.886365643306323*^9, 
   3.886366038222907*^9, 3.886366073647285*^9, 3.886366130018861*^9, 
   3.886366229108131*^9, {3.88636626482434*^9, 3.886366347838192*^9}, 
   3.886366593411806*^9, 3.886368584824757*^9, 3.886368692832632*^9, {
   3.886368816952091*^9, 3.8863688310916986`*^9}, {3.886368937698862*^9, 
   3.886368946734441*^9}, 3.886368980928012*^9, {3.886369661824843*^9, 
   3.88636968643139*^9}, 3.886372266368836*^9},
 CellLabel->
  "Out[545]=",ExpressionUUID->"7ef3bc22-13c9-4d37-b6b4-862f4b7a1360"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "[", "0.55", "]"}]], "Input",
 CellChangeTimes->{{3.886245709970716*^9, 3.88624571748841*^9}, {
  3.886266194904578*^9, 3.886266200614443*^9}},
 CellLabel->
  "In[546]:=",ExpressionUUID->"94a1c09b-e57e-4b4f-be08-8aef27bafd40"],

Cell[BoxData["10.746078283213222`"], "Output",
 CellChangeTimes->{
  3.886245724955971*^9, 3.886245816253265*^9, 3.8862661736647673`*^9, 
   3.8862662060341663`*^9, 3.8862668317151546`*^9, 3.886365643316305*^9, 
   3.886366038234811*^9, 3.88636607367115*^9, 3.886366130031261*^9, 
   3.8863662291184072`*^9, {3.886366264835487*^9, 3.886366347848201*^9}, 
   3.8863665934233923`*^9, 3.8863685848520823`*^9, 3.886368692842907*^9, {
   3.886368816959437*^9, 3.886368831098937*^9}, 3.8863689467427588`*^9, 
   3.886368980936646*^9, {3.8863696618369093`*^9, 3.886369686440624*^9}, 
   3.8863722663764563`*^9},
 CellLabel->
  "Out[546]=",ExpressionUUID->"4f5ea8f2-221f-45c1-a179-c94b8a915e08"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "[", "0.6", "]"}]], "Input",
 CellChangeTimes->{{3.886266808737892*^9, 3.8862668250314302`*^9}, {
  3.8863722337677603`*^9, 3.886372234425582*^9}},
 CellLabel->
  "In[547]:=",ExpressionUUID->"d0d75614-4486-45b8-8e9d-a773bc74cd41"],

Cell[BoxData["10.000000000000043`"], "Output",
 CellChangeTimes->{
  3.8862668317350683`*^9, 3.8863656433325243`*^9, 3.886366038254722*^9, 
   3.8863660736957417`*^9, 3.8863661300544577`*^9, 3.886366229143902*^9, {
   3.886366264856803*^9, 3.886366347868252*^9}, 3.8863665934442463`*^9, 
   3.886368584858449*^9, 3.886368692863385*^9, {3.886368816980013*^9, 
   3.886368831115343*^9}, 3.8863689467635508`*^9, 3.886368980957546*^9, {
   3.886369661863037*^9, 3.88636968646021*^9}, 3.886372266396175*^9},
 CellLabel->
  "Out[547]=",ExpressionUUID->"29142c2b-43cd-445b-8eb3-8530b3df30e6"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "[", "0.65", "]"}]], "Input",
 CellChangeTimes->{
  3.8863687761312113`*^9, 3.886368827067439*^9, {3.88636965308169*^9, 
   3.886369681560944*^9}},
 CellLabel->
  "In[548]:=",ExpressionUUID->"be5955c1-4d77-4078-b321-338a728a2876"],

Cell[BoxData["19.109529749704375`"], "Output",
 CellChangeTimes->{{3.886369661872342*^9, 3.8863696864666443`*^9}, 
   3.8863722664024553`*^9},
 CellLabel->
  "Out[548]=",ExpressionUUID->"031724bf-9e2c-4265-bd1e-778dba265d18"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "[", "0.7", "]"}]], "Input",
 CellChangeTimes->{{3.886372247668548*^9, 3.886372248025392*^9}},
 CellLabel->
  "In[549]:=",ExpressionUUID->"7c26cd6b-050c-4542-8ac5-43d2c20a9fb8"],

Cell[BoxData["74.98942093324564`"], "Output",
 CellChangeTimes->{3.886372266421529*^9},
 CellLabel->
  "Out[549]=",ExpressionUUID->"c0eb188d-cd89-40be-bb56-9fe0a717609f"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "[", "0.75", "]"}]], "Input",
 CellChangeTimes->{{3.8863722569004707`*^9, 3.886372260970395*^9}},
 CellLabel->
  "In[550]:=",ExpressionUUID->"89aeb8b9-36f6-4e66-a541-cfd8d272d3b3"],

Cell[BoxData["604.2963902381334`"], "Output",
 CellChangeTimes->{3.886372266427904*^9},
 CellLabel->
  "Out[550]=",ExpressionUUID->"7f653c95-5318-44a4-9b5f-93fc280403e6"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.8866883667669353`*^9, 
  3.886688368674171*^9}},ExpressionUUID->"c5f5b248-88b6-4e46-9e54-\
bd32eb481465"]
},
WindowSize->{829, 580},
WindowMargins->{{0, Automatic}, {Automatic, 1}},
FrontEndVersion->"12.2 for Mac OS X x86 (64-bit) (December 12, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"c5ced723-74ca-406b-8545-5836840172cb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 325, 5, 30, "Input",ExpressionUUID->"2ba69d38-9fbf-46b8-b9d4-5eaf22c5df35"],
Cell[886, 27, 2401, 38, 1093, "Text",ExpressionUUID->"3f3086f1-6c1a-45b0-9b94-de14f6f460d0"],
Cell[3290, 67, 2164, 44, 52, "Input",ExpressionUUID->"179b6fa4-09df-4050-896c-1e873474453a"],
Cell[5457, 113, 247, 4, 30, "Input",ExpressionUUID->"c4ddf92d-7740-4248-9d94-e0085be039d5"],
Cell[CellGroupData[{
Cell[5729, 121, 359, 9, 30, "Input",ExpressionUUID->"939136d3-589b-4ee0-8a45-d858bfc87826"],
Cell[6091, 132, 4769, 97, 233, "Output",ExpressionUUID->"55a4c815-d58f-4dda-9775-7cb5ab561b3e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10897, 234, 244, 5, 30, "Input",ExpressionUUID->"df530b00-027d-4520-bd6d-95f51302ed02"],
Cell[11144, 241, 1229, 18, 34, "Output",ExpressionUUID->"9c087279-646f-4c5f-b2f6-56516af28d56"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12410, 264, 706, 18, 30, "Input",ExpressionUUID->"7bc97526-b8d6-424c-b9e7-d44e6eff27c2"],
Cell[13119, 284, 7226, 139, 233, "Output",ExpressionUUID->"a66dcbea-304f-4b4e-98b4-ebe2997a1f88"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20382, 428, 603, 15, 30, "Input",ExpressionUUID->"28d0548e-1f8d-405a-9184-5bb1b9d5240d"],
Cell[20988, 445, 6999, 136, 233, "Output",ExpressionUUID->"9f030a14-c0a0-4068-940f-74c6a781cb23"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28024, 586, 603, 15, 30, "Input",ExpressionUUID->"d92a3426-e14d-44e7-be99-e1fa40735cff"],
Cell[28630, 603, 7023, 136, 233, "Output",ExpressionUUID->"b952f162-66be-4bc4-b8c1-82bab1839b72"]
}, Open  ]],
Cell[35668, 742, 248, 6, 35, "Text",ExpressionUUID->"41a96361-fc1b-453f-89b4-e498b0e27f52"],
Cell[CellGroupData[{
Cell[35941, 752, 413, 12, 52, "Input",ExpressionUUID->"26b42d93-c524-4eaa-8b86-9118817303be"],
Cell[36357, 766, 1198, 21, 44, "Message",ExpressionUUID->"fa08f2ba-3efe-4e3c-8b8d-182fb14a497c"]
}, Open  ]],
Cell[37570, 790, 669, 15, 127, "Text",ExpressionUUID->"4870c497-cf11-49b7-9695-f62d4fa430c4"],
Cell[38242, 807, 224, 4, 30, "Input",ExpressionUUID->"6fbb851b-6568-4c96-9532-ba8ea625b693"],
Cell[38469, 813, 1115, 33, 52, "Input",ExpressionUUID->"bccf9706-289c-48df-8278-b45694d34d25"],
Cell[CellGroupData[{
Cell[39609, 850, 1179, 30, 73, "Input",ExpressionUUID->"3dfd632f-5932-4054-979a-5d6182690d97"],
Cell[40791, 882, 15995, 292, 234, "Output",ExpressionUUID->"d2aff15a-362d-41e5-b639-ca9190d07350"]
}, Open  ]],
Cell[56801, 1177, 291, 6, 58, "Text",ExpressionUUID->"09fcf3ab-1b23-4213-8534-9c402e6a4bf5"],
Cell[57095, 1185, 1073, 25, 33, "Input",ExpressionUUID->"b2084bfc-6741-4b3a-a234-020db4a398fa"],
Cell[CellGroupData[{
Cell[58193, 1214, 795, 19, 33, "Input",ExpressionUUID->"38fed2e0-16dd-40b9-84b7-1aa1d91d9f07"],
Cell[58991, 1235, 21677, 372, 223, "Output",ExpressionUUID->"7ef3bc22-13c9-4d37-b6b4-862f4b7a1360"]
}, Open  ]],
Cell[CellGroupData[{
Cell[80705, 1612, 254, 5, 30, "Input",ExpressionUUID->"94a1c09b-e57e-4b4f-be08-8aef27bafd40"],
Cell[80962, 1619, 687, 11, 34, "Output",ExpressionUUID->"4f5ea8f2-221f-45c1-a179-c94b8a915e08"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81686, 1635, 258, 5, 30, "Input",ExpressionUUID->"d0d75614-4486-45b8-8e9d-a773bc74cd41"],
Cell[81947, 1642, 585, 9, 34, "Output",ExpressionUUID->"29142c2b-43cd-445b-8eb3-8530b3df30e6"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82569, 1656, 258, 6, 30, "Input",ExpressionUUID->"be5955c1-4d77-4078-b321-338a728a2876"],
Cell[82830, 1664, 225, 4, 34, "Output",ExpressionUUID->"031724bf-9e2c-4265-bd1e-778dba265d18"]
}, Open  ]],
Cell[CellGroupData[{
Cell[83092, 1673, 205, 4, 30, "Input",ExpressionUUID->"7c26cd6b-050c-4542-8ac5-43d2c20a9fb8"],
Cell[83300, 1679, 170, 3, 34, "Output",ExpressionUUID->"c0eb188d-cd89-40be-bb56-9fe0a717609f"]
}, Open  ]],
Cell[CellGroupData[{
Cell[83507, 1687, 208, 4, 30, "Input",ExpressionUUID->"89aeb8b9-36f6-4e66-a541-cfd8d272d3b3"],
Cell[83718, 1693, 170, 3, 34, "Output",ExpressionUUID->"7f653c95-5318-44a4-9b5f-93fc280403e6"]
}, Open  ]],
Cell[83903, 1699, 154, 3, 30, "Input",ExpressionUUID->"c5f5b248-88b6-4e46-9e54-bd32eb481465"]
}
]
*)

(* End of internal cache information *)

